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> 1 ABSTRACT 



I review recent work, performed in collaboration primarily with N. Read and Jinwu Ye, on the 



O 

CO 

properties of quantum antiferromagnets in two dimensions. The emphasis is on the properties 
of the antiferromagnet in states which do not have any long-range magnetic order. The 
universal spin dynamics in the quantum critical region of number of frustrated and random 
1 I antiferromagnets is studied; indentions for neutron scattering experiments in the Hght.y- 
doped cuprates are noted. The nature of the quantum-disordered phase of non-random 
frustrated antiferromagnets is examined in some detail: the states found have {%) collinear 
spin correlations, spin-Peierls or valence-bond-solid order, and confined spinons, order and 
confined spinons or {%%) coplanar spin correlations, no spin-Peierls order and deconfined 



bosonic spinons. 



1. INTRODUCTION 

The study of quantum antiferromagnets goes back to very early work of Bethe M, Ander- 
son |2| and Kubo [^]. However the subject has recently seen a tremendous amount of renewed 
theoretical and experimental interest. We will restrict the discussion here to antiferromagnets 
described by the following simple Hamiltonian: 

H = Jij&i ■ Sj (1) 

where the S, are spin S quantum spin operators on the sites, i, of a 2-dimensional lattice. 
The antiferromagnetic exchange constants, Jy > 0, are short-ranged and possibly random. 
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The early work [|I], |2], || focused on the nature of the ground states of H with magnetic 
long-range-order (LRO): these are states in which there is a spin-condensate: 

(S<) = m ?; ^ (2) 

The situation changed dramatically with the suggestion by Anderson Q in 1987 that quan- 
tum disordered ground states of two dimensional quantum antiferromagnets are relevant to 
the phenomenon of high temperature superconductivity in the cuprates || . The magnetically 
disordered states are those in which 

(Si) = (3) 

We emphasize that the spin-glass ground state is not of this type; rather it is a magnetically 
ordered ground state in which nij is a random function of i. The following questions arise 
immediately in the study of any quantum disordered phase: 

• Is the quantum-disordered ground state characterized completely by Eqn. (§) ? Or are 
there additional broken symmetries that must be present ? Naively, one may conclude 
that the first of these is the most natural possibility. Many investigators have addressed 
this issue and found however that other broken symmetries appear more often than 
was initially expected. Only under special conditions, which N. Read and I delineated 
in Refs. || 0, ||, does a featureless spin-fluid state appear. This work will be reviewed 
in Section 6. 

• What is the nature of the excitations above the ground state ? Do the individual quanta 
of excitations carry fractional spin quantum numbers (S=l/2) ? Such excitations are 
conventionally called spinons. A second possibility is that the spinons are always 
confined in pairs, leading to excitations which have only integral spin. 

• If the spinons are unconfined, what are their statistics ? In non-random antiferromag- 
nets, the spinons are usually separated from the ground state by a gap, and physical 
conditions with a sufficiently dilute concentration of spinons can always be defined to 
make the question of their statistics meaningful. 

Most of the questions above will be addressed in the course of these lectures. We will use 
a combination of semiclassical and and a number of different large- iV expansions to address 
these issues. In particular, much will be learned by comparing and reconciling the results of 
the different expansions. 

There are several motivations for addressing these and other related questions: 
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The cuprates undergo a transition to a quantum disordered ground state upon doping 
with a small concentration of mobile holes. In La2- x Sr x CuO^ this occurs somewhere 
between x = 0.02 and x = 0.07. At lower dopings the holes are localized and a 
Hamiltonian as in ([I]) is adequate. At larger dopings the holes become mobile, requiring 
study of quantum-disordered antiferromagnets in the presence of propagating holes A 
number of neutron scattering experiments have given us a detailed picture of the spin 



dynamics of these antiferromagnets P, [10], 11, 12, 13, 14, 15, 16], much of which is 
poorly understood. 

The spin-1/2 antiferromagnet on a kagome lattice may posses a quantum disordered 
phase. Recently there have been a number of experimental realizations of Heisenberg 
antiferromagnets on a kagome lattice: the layered kagome antiferromagnets SrCr s _ x - 
Ga i+x 19 , KCr 3 (OH) 6 (S0 4 ) 2 0, HI and layers of 3 He on graphite |(| gl . 



As a noted above, there have been a number of detailed dynamic neutron scattering 
experiments on lightly-doped cuprates ||, [ll], [LI], [12], [13|, [14], [15], [Uj and layered frustrated 



antiferromagnets |20| , |21] , |TT| , [18], [19| . One feature of the dynamic neutron scattering on the 
lightly-doped cuprates (La 2 - x Sr x CuO i for 0.02 < x < 0.05) that is particularly intriguing is 
that the overall frequency scale of the spin excitation spectrum appears to be given simply 
by the absolute temperature. In particular, it appears to be independent of all microscopic 
energy scales e.g. an antiferromagnetic exchange constant. 



We have recently proposed |22] that this anomalous dynamics is a very general property 



of finite T, spin fluctuations in the quantum-critical ]23| region. This is the region in the 



temperature-coupling constant plane where the spin fluctuations are controlled by the crit- 
ical fixed point between the magnetically ordered and the quantum disordered phase. The 
antiferromagnet notices the thermal deviation from T = at a scale shorter than that at 
which it notices that the coupling constants are not exactly critical. The system is neither 
fully ordered or fully disordered at long distances, but remains in limbo controlled by the 
critical point between them; the dynamics exhibits features of both phases and the energy 
scale is set solely by the absolute temperature. In our work [22| universal scaling func- 



tions and exponents for the quantum-critical dynamics of frustrated, doped, and random 
antiferromagnets have been proposed. Complete scaling functions and exponents have been 
calculated for a model system - non-random, frustrated two-dimensional Heisenberg anti- 
ferromagnets with a vector order-parameter. These results will be reviewed in Sections 4 
and 5 and represent the first calculation of dynamic scaling functions in a two-dimensional 
quantum phase transition. 
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Most of the work reviewed here will be on non-random two-dimensional antiferromagnets. 
There will however be a limited discussion of the consequences of randomness in Section 5. 
For the most part, random quantum antiferromagnets are only very poorly understood; there 
has however been some recent progress |24], |25[ for which the reader is invited to consult the 



original papers. 

We will begin this review with a derivation of the coherent-state, path-integral formu- 
lation of quantum antiferromagnets. The representation of the semiclassical fluctuations of 
antiferromagnets by the non-linear sigma model will be discussed in Section 3. In Section 
4, the properties of the O(M) non-linear sigma model will be studied in a large=M theory, 
with a particularly emphasis on the quantum-critical region. This will lead to a more general 
discussion of quantum-criticality in Section 5. Finally the quantum disordered state will be 
discussed in Section 6. 

2. COHERENT STATE PATH INTEGRAL 

In this section we shall present a derivation of a coherent state path integral formulation 
of the dynamics of a generic antiferromagnetic. There will be no approximations made and 
the results are essentially exact. As identical manipulations have to be performed on all the 
spins of the antiferromagnet, we will drop the site index in the following and simply present 
the path integral formulation for a single spin. The path integral for the lattice problem 
will involve direct copies of the single site result. The results of this section are "well- 



known" , although the formulation here was first discussed in Ref. [ 26| ; it has the advantage 
of explicitly preserving spin-rotation invariance. 

As a first step, we need to define the spin coherent states The reader is urged at this 
point to consult the very instructive discussion on the properties and uses of coherent states 
in general in the introduction to the book edited by Klauder [p7] . We restrict our discussion 



here to the case of SU(2) spins. The conventional complete basis for the spin states at each 
site is the usual one corresponding to the eigenstates of S 2 and S z : 

\S, m) with m = — S . . . S (4) 

In the path-integral formulation it is more convenient to use instead an overcomplete basis 
of states |N) labeled by the points N on the surface of the unit sphere. These states are 
normalized, (N|N) = 1 but are not mutually orthogonal. Their crucial properties are that 
the expectation value of the spin is given by 

(N|S|N) = 5N, (5) 
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and the completeness relation 

/^|N)(N| = 1, (6) 

where the integral is over the unit sphere. For N = (0,0,1), the state |N) is easy to 
determine; we have 

|N = (0,0,l)) = |5 , ,m = 5) = |*o) (7) 

We have defined this particular coherent state as a reference state |\&o) as it will be needed 
frequently in the following. For other values of N we can obtain |N) by a SU(2) rotation of 
|^o)- It is n °t difficult to show in this manner that 

|N) =exp(zS + -z*S^) |^o) (8) 

where the relationship between the complex number z and N is simplest in spherical co- 
ordinates: 

N = (sin 6 cos 0, sin 9 sin 0, cos 6) (9) 

Q 

z = --exp(-i0) (10) 

The validity of the Eqn. (|7]) determining the expectation value of S can be verified in a 
straightforward manner. 

It will be useful for our subsequent formulation to rewrite the above results in a somewhat 
different manner, making the SU(2) symmetry more manifest. Define the 2x2 matrix of 
operators S by 

S = \ - ^ z ^ x ~~^ v ) (11) 
\ S x + iSy S z J 

Then Eqn. ([?D can be rewritten as 

(N|4 fl |N) = SQ aP , (12) 

where the matrix Q is 

n ( N * N x -iN y \ _ XT _ , . 

Q= [N a + iN y -N z y )= N -° W 

where a are the Pauli matrices. Furthermore there is a simple relationship between Q and the 
complex number z. In particular, if we use the spin-1/2 version of the operator in Eqn. @ 

f = exp[(_°, J)] (14) 
(U is thus a 2 x 2 matrix), then we find 

Q = Ua z U j (15) 
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We now proceed to the derivation of the coherent state path integral for the partition 
function 

Z = Trexp(-/?ft(S)) (16) 

where we have emphasized that the Hamiltonian ri is a function of the spin operator S; we 
will restrict the following discussion to Hamiltonians in which TC is a linear function of any 
given S on a fixed site. The ri in Eqn. ([!]) is certainly of this type. The transformation of 
Z into a path-integral begins with a standard procedure: we will omit the details of some 



steps, referring the reader to the introductory article in the book edited by Klauder [27 



Briefly, the exponential in Eqn. (|16|) is written as the Trotter product of a large number of 
exponentials each evolving the system over an infinitesimal imaginary-time interval At; the 
identity @ is inserted between all the exponentials and the matrix elements evaluated using 
(§). This yields the result may now be used to obtain the following representation for the 
partition function 



J PN(r)exp jjT dr 



(N(r)\^\N(r))-H(SN(T)) 



(17) 



where ri(SN) is obtained by replacing every occurrence of S in the Hamiltonian by S~N. 
The first term in the action is the Berry phase term, Sb, and represents the overlap between 
the coherent states at two infinitesimally separated times. It can be shown straightforwardly 
from the normalization condition (N|N) = 1 that Sb must be pure imaginary. In the 
remainder of this section we will manipulate Sb into a physically more transparent form 
using the expressions above for the coherent states. 

Clearly, the r-dependence of N(r) implies a r dependent z(r) through fllPl). From (§) 
we have therefore 

i N (r)) = ^-ex P (z(T)S + -z*(r)s) |¥ > (18) 



dr dr 

Taking this derivative is however not so simple: I remind you that if an operator M does 
not commute with its derivative dMj dr then 

4- exp(M) ^ ^ exp(M) (19) 
dr dr 

A careful analysis in fact leads to the general result |28[ 

-f-exp(M) = / 1 rfwexp(M(l - w))^ exp(Mw), (20) 
dr Jo dr 

where u is just a dummy integration variable. Using (|1^) and (|20| ) we find 

Sb = /V<N(t)|A| N (t)> 
Jo dr 

rP r 1 / dz - dz* - \ 

= J Q dr J Q du(N(r,u)\ i—S + - —S\ \N(r,u)) (21) 
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where N(r, u) is defined by 

|N(t,u)) = exp (u (z(r)S + - z*(t)SJ)) |* > (22) 
From this definition, three important properties of N(r, w) are immediately apparent 

N(r,u = l) = N(r) 

N(r,u = 0) = (0,0,1) (23) 
N(r, w) moves with u along the great circle between N(r, u = 0) and N(r, m = 1) 

We can visualize the dependence on u by imagining a string connecting the physical value 
of N(r) = N(t, m = 1) to the North pole, along which u decreases to 0. We can also define 
a u-dependent Q(t,u) from Eqn. fll3"|); furthermore, if we choose 

" ( -°- o )1 < 24 > 



£/ (r, m) = exp 



then the relationship (|15j ) remains valid for all u. Now we use the expression (|12|) for the 
expectation value of S in any coherent state to obtain 



dz dz* 
—Q 2 i(t,u) - -t^Qi 2 {t,u) 



(25) 



As everything is a periodic function of r, we may freely integrate this expression by parts 
and obtain 



S 



B 



Si dr I duTr 





-z*(t) 



z(t) 




d T Q(r,u) 



(26) 



where the trace is over the 2x2 matrix indices. The definitions (ITq ) and (|24|) can be used 
to easily establish the identity 

( Z ( T ) \ 



which when inserted into 



Si 



-z*(r) o ; = -2 g(r ' u) ^^' 

?B) yields the expression for Sb in one of its final forms 



dr I du 



o 



S 



Tr Q(t, u 



dQ(r,u) dQ(T,u)' 



(27) 



(28) 



du dr 

An expression for Sb solely in terms of N(r, u) can be obtained by substituting in (|13"D; this 
yields finally: 



P 

Sb = iS I dr 



du N 



<9N dN 

du dr 



(29) 



This expression has a simple geometric interpretation: it is simply iS times the area on 
the surface of a unit sphere swept out by the string connecting N(r) to the North Pole - 
this is also the oriented area contained within the closed loop defined by N(r), < r < /3, 
N(0) = N(r) on the sphere. Note that the u dependence has dropped out of this result 
which only depends on the values of N(r, u = 1). 
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3. NON-LINEAR SIGMA MODEL 



We shall now obtain a long-wavelength, large-S 1 description of the dynamics of <i-dimensional 
antiferromagnets in the vicinity of a ground state with collinear long-range Neel order. This 
description takes the form of an imaginary-time functional integral over an action which 
includes a non-linear sigma model field theory in d + 1 dimensions; this result was first noted 



by Haldane [29, [30], 31]. The advantage of this approach is that it is the physically most 
transparent, and often the simplest, way of describing many of the results. Furthermore, 
although the derivation depends upon a semiclassical, large-S* limit, all of the results have 



also been obtained in alternative large- iV methods [26, 34, 35] which are valid well into 
the quantum-disordered phase. The chief disadvantages of the present approach are that 
it connect simply describe the appearance of spin-Peierls order in the quantum-disordered 
phase, and it does not have a simple extension to a description of the phases with short-range 
or long-range incommensurate spin correlations 

Let us consider an antiferromagnet on a <i-dimensional hypercubic lattice with only a 
nearest-neighbor exchange interaction J. In the semiclassical limit, the spins will predom- 
inantly orient themselves in opposite directions on the two sublattices. Let n(r, r) be a 
continuum field of unit length which describes the local orientation of this Neel ordering - n 
varies slowly on the scale of a lattice spacing, but the values of n on well separated points can 
be considerably different, leaving open the possibility of a quantum disordered phase with 
no long-range spin order. It will turn out to be necessary to also include a component of the 
spins which is perpendicular to the local orientation of the Neel order - this is described by 
the continuum field L(r, r). We have therefore on the site i of the lattice 



N(i) « eXfity 1 - « M L 2 (r,) + a d L( ri ), (30) 

where £j equals ±1 on the two sublattices, a is the lattice spacing, and 

n 2 = l n-L = L 2 <a~ M (31) 

These relationship holds for all values of r and u, and n, L satisfy the same boundary 
conditions in u as N. We insert the decomposition for N into 7Y(N(r)) and expand the 
result in gradients and powers of L. This yields 

H = \J d d r [jS 2 a 2 - d (V r nf + 2dJS 2 a d V 

I I d d v k(V r n) 2 + x±5 2 L 2 l (32) 



2 

the second equation defines the thermodynamic spin- wave stiffness, p s , and the transverse 
susceptibility x±- If we na d used a different form for TC with modified short-range exchange 
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interactions, the continuum limit of Ti would have been the same but with new values of p s 
and x±- 

To complete the expression for the coherent state path-integral in the the continuum 



limit we also need the expression for S B in terms of n, L. We insert (30) into the ( p9| ) and 
retain terms upto linear order in L: this yields 



S B = S' B + iS [ d d r ( dr f 1 du 
J Jo Jo 



n ■ 



' dn dL 

du dr 



n 



' dL dn 

du dr 



dn dn 

du dr 



where 



S', 



B 



P 



dr \ du n, 



<9n, : diij 

x 



(33) 



(34) 



du dr 

The continuum limit of S' B is not simple due to the presence of the rapidly oscillating 
prefactor 6{. In two spatial dimensions, a careful examination by several investigators p2| 
showed that this term is identically zero for all smooth spin configurations. However when 



singular spin configurations are permitted, there can be a net Berry phase |p3| ; this will 
studied in Section 6.C. We also note in passing that in one-dimensional antiferromagnets p9[ 
it is S' B which gives rise to the topological (9-term. 

The remainder of Sb can be simplified further. Using the fact that the vectors L, dn/dr, 
dn/du are all perpendicular to n and hence lie in a plane and have a vanishing triple product 
we find 



S B = S' B + iS f d d r [ dr C du 
J Jo Jo 



d 



dn 



d 



n 



L *9 



(35) 



dr V \du 

The total r derivative yields as all fields are periodic in r, while the total u derivative 
yields a surface contribution at u — 1. This gives finally 



S, 



S' B -iS I d d r 



drL • n x 



dn" 



(36) 



Putting together ( |32|) and (|36|) we obtain the following continuum limit path-integral for the 
ant if erromagnet 

VnVL exp(5' n ) 



Z 
S n 



S', 



B 



1 

2 Jo 



P 



dr / d d r 



Ps (V r n) 2 + X ±S 2 L 2 + 2iSL ■ L x ^ 



(37) 



The functional integral over L is simply a gaussian and can therefore be carried out 



(V r n) 2 + ~{d T n) 



(38) 



Sn = S ' B + \So dT I 

where the spin- wave velocity c = ■Jp~sXX- This is the action for ad+1 dimensional non-linear 
sigma model with a residual Berry phase term. As already noted, the Berry phase terms 
only makes a non-zero contribution for topologically non-trivial spin configurations. 
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4. O(M) NON-LINEAR SIGMA MODEL FOR LARGE M 



We will now study the properties of the O(M) non-linear sigma model field theory in a 1/M 
expansion; such an expansion was carried out for the two-dimensional model by Polyakov ||36|| . 
The results in this section represent work carried out with Jinwu Ye |22, |37|j . 



As we have seen in the previous section, the 0(3) non-linear sigma model, combined with 
some additional Berry phase terms, describes the long-wavelength dynamics of the SU(2) 
antiferromagnet in the semiclassical limit. The Berry phases will be completely ignored in 
this section - their consequences will be considered later. Furthermore, the 0(M) model 
for M > 3 is not simply related to any quantum antiferromagnet. Our main motivation 
for looking at this model is that, with proper interpretation, it presents the simplest way at 
arriving at some of the basic results and exploring crude features of the phase diagram of the 
SU(2) antiferromagnet. The omited Berry phases do induce important differences between 
the SU(2) antiferromagnet and the 0(3) non- linear sigma model - these will be explored in 
subsequent subsections. 

In this section, we will deal exclusively with the following non-linear sigma model field 
theory 

Z = Jvn a S(nl-l)exp(~ J d 2 r df [(V r n a ) 2 + (d f n a ) 2 ] J , (39) 



where the index a runs from 1 to M. This action can be obtained from the semiclassical 
action for the antiferromagnet ([53) by omiting the Berry phase term S' B , and introducing 
the rescaled time f = cr. We will henceforth omit the tilde on the r and use units in which 
c = 1 - it is trivial to reinsert appropriate factors of c in the final results. The coupling 
constant g is given by 

Mc 
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Ps 

= — for the large S SU(2) antiferromagnet in d — 2 (40) 
The large M analysis of Z begins with the introduction of the rescaled field 

h a = VMn a (41) 
and the imposition of the constraint by a Lagrange multiplier A. This transforms Z into 

Z = J Vn a VX exp J d 2 r ^ dr [(V r n a ) 2 + (d T n a ) 2 + i\(n 2 a - M)fj , (42) 

This action is quadratic in the n a , which can therefore be integrated out. This induces an 
effective action for the A field which has the useful feature of having all its M dependence in a 



10 



prefactor. Therefore, for large M the A functional integral can be evaluated by an expansion 
about its saddle point. 

Let us look more carefully at the solution at M = oo when the saddle point action is the 
exact answer. At the saddle point we parametrize 

= m 2 (43) 

where m is a function of g, T, and an upper cutoff which will be determined below. The 
correlator of the n a field is the stagered dynamic spin susceptibility and is given by 

g 

X(k, u n ) = (n a (k, u n )n a (-k, -u n )) = — - — — (44) 

k z + aj£ + m z 

where we have introduced the momentum k and the Matsubara frequency u n . Neutron 
scattering experiments yield a direct measurement of the imaginary part of x {x") f° r rea l 
frequencies. At M = oo, analytic continuation of the above result yields 

X"(k, u) = —, = | = = (S(u - Vk 2 + m 2 ) - 5{uj + Vk 2 + m 2 )) (45) 
2vk 2 + m 2 v ' 

Another quantity often acessed in neutron scattering is the momentum integrated local 
susceptibility, xl, which is given by 

X ' L{U) = I 4S X " (k ' U) = 4 {6{UJ ~ m) ~ d{ ~ U " m)) (46) 
where 6{x) is the stop function which is non-zero only for positive x. An important short- 
coming of the M = oo result is the absence of any damping which is always present at finite 
T - this only appears at order 1/M which will be discussed later. The damping will broaden 
the delta function peaks in x(k, u) and fill in the gap region (\cu\ < m) in Xl(^) at all 
finite temperatures and all values of g - true delta functions and gaps can however survive at 
T = 0. We can also examine the spin structure factor, S(k) which is the equal-time spin-spin 
correlation function in momentum space 

S(k) = T^2(n a (k,u n )n a (-k,-uj n )) 



2Vk 2 + m? V 2T 
The spin-correlation length, £, is therefore given by 

.2 1 dS 



1 coth ( Vk2 t m2 ) (47) 



k=0 



S(0) dk 2 
1 (\ m/T \ 

^2 \ 2 + 2 sinh(m/T) I ^ 
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The factor multiplying l/m 2 decreases monotonically from 1 to 1/2 as m/T increases from 0. 
Thus m is essentially the inverse spin correlation length, apart from an innocuous numerical 
factor between 1 and l/\/2. 

We now determine m. The saddle-point equation is simply the constraint n 2 = 1, or 
more explicitly 

TW— 9 ~ = 1 (49) 

i^J 4tt 2 k 2 + cu 2 + m 2 1 ' 

where T — 1/(5 is the temperature (we are using units in which fc B — fr — 1). It is easy to see 
that this equation is divergent in the ultraviolet, and it is therefore necessary to introduce 
an ultraviolet cut-off in the momentum integration. The most convenient method is to use 
a Pauli-Villars cut-off A. This transforms the constraint equation into 

d 2 k f g g 



T 



4tt 2 (k 2 + cu 2 +m 2 k 2 + cu 2 +A 2 ) 1 (50) 



The momentum integration can be carried out and yields 



4vr \u*+m?J 
The frequency summation can be done exactly to give 



gT / sinh(A/2T) \ _ 

27r ln l v sinh(m/2T)J _i (52j 
Finally, we can solve for the dependence of m on T, g, and A 

m = 2TArcsinh ^exp ^ — sinh (^Tj^J^j (53) 

A remarkable amount of information is contained in this innocuous looking equation. By 
examining the T — > limit of this equation, it is immediately apparent that the behavior of 
m is quite different depending upon whether g is smaller, larger, or close to a critical value 
g c given by 

An 

9c = T (54) 
We examine the three different phases separately. 

Ordered Phase - g < g c 

For g < g c we find for small T that 

m ~ r 1 ~ 2Texp (-^ - -)) (55) 



.9 9 
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Thus the correlation length diverges exponentially as T — > and long-range-Neel order 
appears at T = 0. A renormalization group analysis of the non-linear sigma model by 



Chakraborty et. al. [23| yielded the same functional dependence of the correlation length on 
T; this agreement gives us some confidence on the utility of the present large M expansion. 
We refer the reader to Ref. [23] for further discussion on the dynamic properties of the 
ordered phase - the semiclassical approach discussed there is a little more convenient for 
analyzing this region. 

Quantum Disordered Phase - g > g c 

We now find for small T that 

1 Aix(g - g c ) 

m ~ — ~ 56 

2£ gg c K } 

with exponentially small corrections at low T. In this case there is a finite correlation length 
at T = 0. Furthermore the excitation spectrum has an energy gap = cm. The low-lying 
excitations are the massive spin-1 (for 0(3)) n a quanta -they are thus triply degenerate. 
There will also be a spinless collective mode with a gap represented by the fluctuations of A 
about its saddle-point value. 

As g approaches g c we expect on general scaling grounds that 

Z~\g-9c\- v (57) 



Comparing with (p6|) we have v — 1 at M — oo. We have considered the 1/M corrections to 



this result. The structure of the perturbation theory is very similar to the 1/M expansion 



for soft-spin Ginzburg Landau models as discussed in the book by Ma |38[] ; we will therefore 
omit the details. We find 

v=\- (58) 

Critical Region - g « g c 

In this region both Sg = g — g c and T are small, and their ratio can be arbitrary. Before 
examining the result ( |53"D let us see what can be predicted on general scaling grounds for 
general M. At g — g c and T = we have a quantum field theory at its critical point. 
Turning on a finite T is equivalent to placing this critical system in a box which is finite 
in the imaginary time direction, while tuning g away from g c at T = induces a finite 
correlation length which scales like (|5T|). The behavior of £ for both T finite and g ^ g c is 
therefore predicted by the principles of finite-size scaling (after restoring factors of h, c, and 
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Figure 1: The universal function X(x) for the critical behavior of the spin-correlation length 
at M = oo. X(0) is close to, but not exactly, unity 



h) 



k B T 



X ai 



sgn(ff ~ 9c)\9 ~ 9c 
T 



V 



(59) 



where a\ is the only non- universal, cutoff dependent quantity and X{x) is a universal func- 
tion. X(x) is a smooth function of x except at x = where it is continuous but not 



in this form with v — 1, a\ = 27r/<^ and X(x) is defined by the following equations 



A plot of the function X(x) is shown in Fig. |l|. 

The quantum-critical region is defined by the inequality \5g\ u /T < 1 (Fig |2|). In this 
region the antiferromagnet notices that it is finite in the time direction at a scale which 
is shorter than which it notices that the coupling g is not exactly g c . Thus the critical 
spin fluctuations are quenched at an energy scale which is determined completely by the 



differentiable. It is not difficult to see that the M = oo Eqns. ( |48| ) and fl53|) can be written 




(60) 
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Figure 2: Phase diagram of TC in two dimensions. The magnetic LRO can he either spin- 
glass or Neel, and is present only at T = 0. The boundaries of the quantum-critical region 
are T ~ \g — g c \ zu . For non-random 7i which have commensurate, collinear, Neel LRO for 
g < g c , all of the quantum- disordered region (g > g c ) has spin-Peierls order at T = 0-this 
order extends to part of the quantum disordered region at finite T. 

temperature and not any underlying antiferromagnetic exchange constant. In this regime 
we may simply put x = in the formulae in the previous paragraph. 

The full dynamic susceptibility also satisfies scaling functions in the critical region. For 
simplicity we will only consider them in the quantum critical region with Sg = 0: the 
extension to finite Sg is straightforward, at least at M = oo. Application of finite size- 
scaling to the dynamic susceptibility yields 

where y) is a universal scaling function of both arguments, 02 is the only non- universal 
quantity, and rj is the usual critical exponent determining the decay of spin correlations at 
criticality. A normalization condition is necessary to fix the scale of $: the most convenient 
is to use 

-a* = 1 ^ 

ux x=0,j/=0 

The scaling form for the local susceptibility xl is a little more subtle. If we perform the 
momentum integration of ([H]) we notice immediately that the integral is ultraviolet divergent 
because for large k, x ~ k _2+,? and 77 > (see below) for the present model. However this 
divergence is present only in the real part of x- The imaginary part of % will only involve 
excitations only on-shell and will therefore have rapidly vanishing spectral weight as |k| 
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becomes much larger than uj/c, leading to a convergent momentum space integral. Therefore 
we have the scaling form 

where F is a universal scaling function and 03 is a non- universal constant. A remarkable 
feature of this result is that the energy scale for all the excitations in the system is set solely 
by the absolute temperature and is independent of any microscopic energy scale. 

The results for the quantum-critical scaling functions at M = 00 can be easily determined 
from the results of this section. We find rj = 0, 

${x,y) = 1 (64) 

x 2 + B 2 — (y + it) 1 

where e is a positive infinitesimal, and 

f(v) = \ (Kv - e) - o(- y - e)) . (65) 

The quantity G is a pure number given by 

6 = /(0) = 2 In (^Y^j (66) 

We have computed these scaling functions to order 1/M f22 |. This is a fairly non-trivial 
calculation and details will not be presented here. We present the results below which 
required about 40 hours of computation on a vectorized computer. Our results for Im<3> and 
F for M = 3 are summarized in Figs. ^| and (|. Note that Im$ had a well defined (at least 
for large |k| ) spin- wave peak at a frequency close to lu — c|k|. The peak is broadened due 
to a universal damping arising from spin-wave interactions. Analytic forms for $ can be 
obtained in various regimes. We have 

ReS)" 1 = C Q 2 + x 2 + ... x,y small (67) 



Co 1 = © ( 1 + ^) (68) 



The universal number Cq, to order 1/M, is 

N 

Im$ has a singular behavior for x, y small: 

1 / 36 2 \ 

Im$(x = 0, y) ~ — exp ( - — 1 (69) 

while 

1 / 30 2 \ . 

Im$(y < x) ~ — yexp I -— ^ I (70) 
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Figure 3: The imaginary part of the universal susceptibility in the quantum-critical region, 
$, as a function of x = hcq/ (k B T) and y = hu/(k B T) for a non-random square lattice AFM 
which undergoes a T = transition from Neel LRO to a quantum- disordered phase. The 
results have been computed in a 1/M expansion to order 1/M and evaluated for M = 3. 
The two-loop diagrams were analytically continued to real frequencies and the integrals 
then evaluated numerically. The shoulder on the peaks is due to a threshold towards three 
spin-wave decay. 
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Figure 4: The imaginary part of the universal local susceptibility, F, for the same model as 
in the previous figure. We have F(y) = y~^ / dx Im§(x,y). The oscillations at large y are 
due to a finite step-size in the momentum integrations. 
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These peculiar singularities are probably artifacts of the large M expansion and occur only 
in the region x, y < l/M - the naive expectation that Im$ ~ y for small y is probably 
correct. With either x, y large, $ has the form 

Do „ 0.3426 

* = F^ + - ; Dq ^ 1 -— (71) 

The exponent r\ has the known p9| expansion 



8 512 

^ = 3^-27^ >0 (72) 
The scaling function for the local susceptibility, F(y), has the limiting forms 

F{y) = sgn(y)°-^\y\^ y « 1 ; F(y) = sm (y) ^^ y»l (73) 

As 7] is small, F is almost linear at small y. Note that the finite temperature l/M fluctuations 
have filled in the gap in the spectral functions at all frequencies. These scaling functions 
represent the universal dissipative dynamics of critical spin wave fluctuations. 

5. QUANTUM CRITICAL DYNAMICS OF 2D ANTIFERROMAGNETS 

The universal dynamic scaling forms discussed in the last section in the context of unfrus- 
trated antiferromagnets can in fact be generalized to describe the quantum critical region 
of a large number of frustrated and random quantum antiferromagnets. Of course, explicit 
calculations of the exponents and scaling functions is much more difficult and only a few 
results are currently available. There has however been some very recent progress [^3, ^5 
which will not be reviewed here. 

Our motivation for considering quantum critical fluctuations in random, frustrated anti- 
ferromagnets is provided by recent neutron scattering experiments on lightly-doped La^- X - 
Sr x CuO^ where 0.02 < x < 0.07. In this doping range it has been noted by the experimen- 
talists 0, |TTJ that the dynamic spin susceptibility measured as a function of frequency, u, 



and temperature T follows the following scaling form with reasonable accuracy 



4 ^'<(k,.) = /( M )^-j (74) 

Experimental results for the functions F, I from Ref. ]TB] are shown in Figs [5] and ^. We 
will argue below that the theory predicts that the prefactor / must be of the form 

/ ~ (75) 
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Figure 5: Neutron scattering results for the scaling function F from Ref. [15] 



Figure 6: Neutron scattering results for the normalization factor I from Ref. [15] 
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where the exponent // is expected to satisy — 1 < fi < 0. A fit of this form by B. Keimer [20 
to the experimental result yielded \x = —0.41 ± 0.05. 

A remarkable feature of the above result for x is that the frequency scale for the spin fluc- 
tuations is set completely by the absolute temperature. The underlying exchange constants 



appear to modify only the prefactor of the scaling form. We have recently proposed ^2 



that this scaling of the frequency scale with temperature is a rather general property of spin 
fluctuations of two-dimensional antiferromagnets in their quantum critical region. Let us 
discuss this in the context of the general Hamiltonian Ti which we will use to model the 
lightly-doped cuprates. In the lightly doped region the holes are localized at low tempera- 
tures, indicating that a suitable form of Ti will be sufficient to describe the spins. A specific 



form has in fact been proposed by Gooding and Mailhot |^T| an d yielded reasonable results 
on the doping dependence of the zero temperature spin correlation length. 

We will find it necessary to distinguish between two different types of magnetic LRO: 

(A) Neel LRO, in which case 

m, ~ e iQR * (76) 
with Q the Neel ordering wavevector, and 

(B) spin-glass LRO, in which case m, can have an arbitrary dependence on i, specific to 
the particular realization of the randomness. The lower critical dimension of the Heisenberg 
spin-glass ||2| may be larger than 3 - in this case the spin-glass LRO will not survive to 
any finite T, even in the presence of a coupling between the layers. This, however, does not 
preclude the existence of spin-glass LRO at T = 0. 

Consider now a T = phase transition between the magnetic LRO and the quantum- 
disorder, induced by varying a coupling constant g which is dependent on the ratio's of the 
Jij in Ti. Let the transition occur at a critical value g = g c . If the transition is second-order 
then there will be a diverging scale 

Cg^\9-9c\- v (77) 

which is the distance at which the antiferromagnet first notices that its properties are not 
critical; at distances shorter than £ 9 the spin fluctuations are indistinguishable from those at 
criticality. For Neel LRO, £ 9 is identical at T = to the spin correlation length £ considered 
in Section 4. Therefore l/£ is the width of the peak in the spin structure-factor at the 
ordering wavevector Q. The meaning of £ 9 is more subtle for the case of spin-glass LRO. 
As the condensate rrij is a random function of i, there will be no narrowing of the peak in 
the structure-factor. However it is possible to define a diverging length scale associated with 
correlations the Edwards-Anderson order-parameter qEA', £g will therefore be related to the 



decay of certain four-spin correlation functions [42] 
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At finite temperature and near g = g c we can define a second length scale £t associated 
with the distances at which thermal effects become apparent. The dependence of £r on T 
can be deduced from finite-size scaling. At g = g c the antiferromagnet is a d + 1 dimensional 
critical system which is described by a scale-invariant field theory. However, it is not neces- 
sary in general for the theory to be Lorentz invariant, as distances along the imaginary time 
direction can scale with a non-trivial exponent, z, with respect to spatial distances. The 
exponent z is also referred to as the dynamic critical exponent. As 1/T acts as a finite-size 
in the time direction we therefore expect 

fr ~ rm < 78 > 

Now imagine that the system is in a region where £t < £,g (we assume implicitly that both 
these scales are much larger than the lattice spacing) which defines the quantum-critical 
region - see Fig. Q The boundaries of this region are therefore specified by 

T~\g-g c r (79) 

Under these circumstances the antiferromagnet will notice the finite temperature before it 
notices the deviation from the purely critical behavior at g — g c . The spin-fluctuations will 
therefore be controlled by the repulsive flow away from the critical fixed point, cutoff by a 
finite scale which is set by the temperature. Both these properties are universal features of 
the critical fixed point and the dynamic spin susceptibility is therefore expected to obey a 
universal scaling form. We will study the scaling form separately for the transition from the 
two different types of magnetic long-range-order distinguished above: 

5. A Neel order 

This section generalizes the results of Section 4 to include the effects of disorder. We em- 
phasize however that the disorder is not strong enough to destroy the small-g Neel ordered 
phase. 

The first question which must addressed is whether an infinitesimal amount of disorder 
is relevant at the fixed point describing the transition in the pure system. We begin by 
assuming that disorder is irrelevant. In this case the critical exponents will be unmodified. 
In particular the exponent v will preserve the value v = 0.705±0.005 of the pure system. 



However by a result of Chayes et.al. H] the value of v in any random system must satisfy 

v > ^ (80) 

where d is the dimensionality of the space over which disorder is uncorrelated. In our case, 
disorder is uncorrelated only along the two spatial directions, which leads to v > 1. As 
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this inequality is violated by the above value of v, our hypothesis is inconsistent. Therefore 
disorder is a relevant perturbation and must modify the universality class of the transition. 

The random fixed-point is quite difficult to access in general, but a few of its properties 
can be delineated by general considerations. Firstly, the space-time anisotropy induced by 
disorder implies that the critical-theory need not be Lorentz invariant (unlike Section 4) and 
therefore 

z ± 1 (81) 

The results in Section 4 for the dynamic scaling form for x(k, u) can now be easily gen- 
eralized. At T = the static spin susceptibility, x> wm have a divergence at g = g c and 
wavevector k = Q: 

x(k = Q,^ = 0)~ . 1 . (82) 

\g -9cV 

with 7 = (2 — rj)p. At finite k, u,T finite-size scaling yields the scaling form |22j 

. ai , /aUk — Q| %uj \ 

^ ") = J^T® { T y z . k^f j (83) 

where ai,a2 are non- universal constants, and $ is a universal, complex function of both 
arguments. The deviations from quantum-criticality lead to an additional dependence of $ 
on £t/£ : this number is small in the QC region and has been set to 0. The scaling form 
for the local susceptibility can be obtained by integrating ( [S3|) over k. The integral over k 
will be convergent in the ultraviolet provided i] < - available results, which are discussed 
below, strongly indicate that all random systems in fact have rj < 0. We may therefore freely 
integrate over k and obtain 

X "M = a 3 M" F (84) 

with 

H = r]/z, (85) 

F(y) = y-> 1 |01m$(f,y) (86) 

a universal function, and 03 a non-universal number. The real part x'l w ^ 8i ^ so obey an 
identical scaling form as long as 77 < 0. The asymptotic forms for F for small and large 
arguments can be deduced from general considerations. For large y, or Tko ^> /c#T, we 
expect x to become independent of T; therefore 

F(y) ~ sgn(y) for \y\ > 1 (87) 

For small u, but T finite, we expect that x'l ~ 00 ~ no anomalous scaling with uo is expected 
at finite temperature. Therefore 

F{y) ~ lyl 1 -^^) for \y\ < 1 (88) 
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We emphasize again that all non-universal energy scales only appear in the prefactor a% and 
the frequency scale in F is determined solely by T. 

Let us now examine two simple realizations of randomness in 7i for which exponents can 
be estimated from the literature. 
(a) Weak bond randomness 

Take the simplest pure model which has a transition from Neel LRO to quantum disorder - 
we will argue in Section 6 that this is the square lattice spin-1/2 antiferromagnet with first 
(Ji) and second (J 2 ) neighbor interactions - the Ji — J2 model. The transition is expected 
to be described by the non-linear sigma model of Section 4. (As noted in Section 6.C, the 
quantum-disordered phase of this model possesses spin-Peierls LRO - it has been argued 



however that the coupling to the spin-Peierls order is dangerously irrelevant |22|, [46 
and does not modify two-spin correlation functions in the quantum-critical region; the spin- 
Peierls order will therefore be neglected here.) Now add a small fluctuation in the J\ bonds 
of the J\ — J2 model : i.e. J\ — + J\ + 5J\ where 5 J\ is random, with zero mean, and r.m.s. 
variance <C J\. This last condition ensures that a Neel LRO to quantum disorder transition 
will continue to occur as a function of J2/ J\. Moreover, as the disorder is weak, the mapping 
to the non-linear sigma model of Section 3 will continue to work, leading now to a modified 
partition function (|39D with the coupling g a random function of r. A soft-spin version 



of this action has been examined by Dorogovstev and Boyanovsky and Cardy [47]. They 
considered a field theory with a M component order parameter which had random couplings 
in d = 4 — e — e T space dimensions and e r time dimensions. The theory was then examined 
in a double expansion in e, e T |47]]. The expansion is poorly-behaved, and for the case of 
interest here (M = 3, e = 1, e T = l) the random fixed-point has the exponent estimates 

77 = -0.17, 2 = 1.21, v — 0.64, //=-0.15. (89) 

Note (i) fi, 7] are negative, unlike the pure fixed point and (ii) v is smaller than 2/d, violating 



the required bound [44]. This latter discrepancy is not a cause for great concern as it is clear 
from the low-order results that the series is badly behaved and that there are large higher- 
order corrections. 

(b) Static holes on square lattice vertices 

Consider next a TC on the square lattice with only J\ couplings, but with a small concentration 
of static, spinless holes on the vertices; this model will display a Neel LRO to quantum- 
disorder transition at a critical concentration of holes. In the coherent-state path-integral 
formulation of the pure model, discussed in Section 3, each spin contributes a Berry phase 
which is almost completely canceled in the continuum limit between the contributions of 
the two sublattices. The model with holes will have large regions with unequal numbers of 
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spins on the two sublattices: such regions will contribute a Berry phase which will almost 
certainly be relevant at long distances. Therefore the field theory of Ref. H7| is not expected 



to describe the Neel LRO to quantum-disorder transition in this case. It appears instead 
to be a model which possesses relevant spin configurations whic have complex weights in 
Euclidean time, and not amenable to a field-theoretic analysis by existing methods. The 
only available results are those of Wan et. al. E5[ who performed a cluster expansion in the 



concentration of spins i.e. they expanded about configurations with the maximum number 
of holes. Their series analysis yielded the exponents 

77 = -0.6, z = 1.7, v = 0.8, n = -0.35. (90) 

Note again that /j, T] < 0, although the violation of v > 2/d suggests problems with the series 
extrapolations. 

5.B Spin-glass order 

We now turn to the case where the degree of frustration is large enough to destroy the Neel 
LRO state and induce first a state with spin-glass LRO. The transition to quantum-disorder 
occurs subsequently from the spin-glass phase. We will argue below that the lightly doped 
cuprates are better described by this scenario. 

Consider then the scaling of x m the quantum-critical region of the spin-glass LRO to 
quantum-disorder transition. As the condensate rrij is a random function of i, we do not 
expect any singular behavior as a function of k. In particular, the scaling form (|83|) will not 
be obeyed. However, the local susceptibility xr,(a> n ) will be quite sensitive to the transition. 
In particular, xl measures on-site spin-correlation functions 

rV(k B T) 

XlM = / dre^Cir ), 

Jo 



C(r) = (8,(0) • S 4 (r)) (91) 

(where the bar represents an average over all the sites i) which acquire long-range order in 
time in the spin-glass phase. Thus in the T = spin-glass phase 



lim CM = mf = q EA > (92) 

r— >oo 

where q E A is the Edwards-Anderson order-parameter. At the phase transition this order- 
parameter will vanish as 

q EA ~{g c -gY (93) 

At T = 0, the dynamic-scaling hypothesis then implies that near g = g c 

C(r) = (g c -gfh(r\g-g c n (94) 
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where h is a universal scaling function which tends to a constant for large argument. For 
T ^ \g — 9c\~ zl ' i C(t) should be independent of g — g c \ a standard argument then implies 
that at criticality 

Using (^T|) it follows then that x'L obeys a scaling form identical to fl84| ) with the exponent 
fi now given by 

= -1 + A (96) 

zv 

The function F continues to satisfy the asymptotic forms in fl57| ) and fl88|). We emphasize 
that the above scaling forms are valid, despite the failure of the k dependent scaling (|83|). 
As everything is expected to be non-singular as a function of k, we in fact expect (|4|) to 
hold independently at all k . 

Until recently |23|, |2"5| , there were essentially no theoretical results for /i for a transition 
in this class. The identity ( |96| ) clearly shows (j, > —1. It is also plausible that /i < 0. 
This follows from numerical studies of the quantum-disordered phase of random, spin-1/2 
antiferromagnets j|9| which found a divergent zero-temperature local susceptibility in the 
quantum-disordered phase. The susceptibility at the critical point is therefore also expected 
to be divergent, implying n < 0. 

Finally, we return to a discussion of the experimental situation. An important property 
of the experiments is that the spin correlation length becomes independent of temperature, 
over a significant range of low temperatures. This is clearly incompatible with the presence of 
a T = transition from Neel LRO to quantum-disorder which would have £ ~ 1/T at g = g c . 
This leaves the remaining possibility of a spin-glass LRO to quantum-disorder transition. 

Further support for this scenario is provided by the work of Gooding and Mailhot [f4T 



they considered a model of the doped cuprates in which the dopant holes are localized on the 
oxygen sites and the Cu—O—Cu complex on that bond has been replaced by a single spin-1/2 
doublet. This model leads naturally to strong frustration in the ground state. A simulation 
of the classical version led a T = correlation length which was in good agreement with the 
experimental results. Further, their results seem to indicate that the in the 2D system, Neel 
LRO is destroyed and spin-glass LRO appears at any non-zero concentration of dopants. 

The main experimental measurements which then need to be understood are the value of 
[i = —0.41 ± 0.05 and the functional form of F. The fact that fi < verifies experimentally 
the theoretical prediction that disorder will destroy the pure fixed point which had fi > 0. 
More detailed theoretical studies of quantum spin glass are clearly called for |24], |25 . 
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6. QUANTUM DISORDERED PHASES OF NON-RANDOM 2D ANTIFER- 
ROMAGNETS 



We now turn to an examination of the structure of the possible quantum-disordered phases 
of 7i in the absence of randomness. Our study of these phases has so far been limited to 
the large M, O(M) non-linear sigma model in Section 4. The structure of the quantum- 
disordered phase found there was rather simple: the ground state was a featureless spin-fluid 
with massive integer spin excitations. No half-integer spin excitations were found. We shall 
find upon careful examination here that this simple picture is almost never valid ||, [7j . The 
only exception will be even-integer antiferromagnets on the square lattice which have a two- 



sublattice Neel ground state in the classical limit ||50|| . We note however that the results on 
the quantum-critical dynamic scaling functions found above are not modified by any of the 
effects to be discussed in this section: the properties of the quantum-disordered phase are 
controlled by a strong-coupling fixed point while quantum criticality is determined by the 
critical fixed point |22| . 

A careful examination of the quantum-disordered phase requires an approach which is 
designed explicitly to be valid in a region well separated from Neel LRO; the non-linear sigma 
model is essentially a semiclassical expansion and can only approach the disordered phase 
from the the Neel state. To this end, we introduce the Schwinger boson description [[51]]. For 
the group SU(2) the complete set of (25* + 1) states @ on site i are represented as follows 

\S,m) = 1 (bl^birnO). (97) 

J (S + my.(S — my. 

We have introduced two flavors of bosons on each site, created by the canonical operator b ia , 
with a =|, |, and |0) is the vacuum with no bosons. The total number of bosons, n b is the 
same for all the states; therefore 

bib? = n b (98) 

with n b = 2S (we will henceforth assume an implied summation over repeated upper and 
lower indices). It is not difficult to see that the above representation of the states is com- 
pletely equivalent to the following operator identity between the spin and boson operators 



Sia = lb\ a a^ (99) 



2 

where a = x,y,z and the a a are the usual 2x2 Pauli matrices. The spin-states on two sites 
i, j can combine to form a singlet in a unique manner - the wavefunction of the singlet state 
is particularly simple in the boson formulation: 

, 25 



O a/3 *4A) io> (ioo) 
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Finally we note that, using the constraint fl98|), the following Fierz-type identity can be 
established 

& s blbf} = -2S, • S j + n 2 b /2 + 5 {j n b (101) 
where e is the totally antisymmetric 2x2 tensor 



1 
-1 



(102) 



This implies that TC can be rewritten in the form (apart from an additive constant) 



This form makes it clear that 7i counts the number of singlet bonds. 

We have so far defined a one-parameter (n&) family of models TC for a fixed realization of 
the Jij. Increasing n& makes the system more classical and a large rib expansion is therefore 
not suitable for studying the quantum-disordered phase. For this reason we introduce a 
second parameter - the flavor index a on the bosons is allowed to run from 1 . . . 2N with N 
an arbitrary integer. This therefore allows the bosons to transform under SU(2N) rotations. 
However the SU(2N) symmetry turns out to be too large. We want to impose the additional 
restriction that the spins on a pair of sites be able to combine to form a singlet state, thus 
generalizing the valence-bond structure of SU(2) - this valence-bond formation is clearly a 
crucial feature determining the structure of the quantum disordered phase. It is well-known 
that this is impossible for SU (2N) for N > 1 - there is no generalization of the second-rank, 
antisymmetric, invariant tensor e to general SU(2N). 

The proper generalization turns out to be to the group Sp(N) 0. This group is defined 
by the set of 2N x 2N unitary matrices U such that 



U T JU = J 



(104) 



where 



•Ja/3 — J 



a/3 



(105) 



is the generalization of the e tensor to N > 1. It is clear that Sp(N) C SU(2N) for N > 1, 
while Sp(l) = SU(2). The bf bosons transform as the fundamental representation of Sp(N); 
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Figure 7: Phase diagram of the 2D Sp(N) antiferromagnet Ti as a function of the "spin" rib 

the "spins" on the lattice therefore belong to the symmetric product of rib fundamentals, 
which is also an irreducible representation. Valence bonds 



^ u ia u ia 



(106) 



can be formed between any two sites; this operator is a singlet under Sp(N) because of ( |104| ). 
The form (|103|) of 7i has a natural generalization to general Sp(N): 



J, 



(107) 



where the indices 7,5 now run over 1 . . . 2N. We recall also that the constraint ( p8|) 
must be imposed on every site of the lattice. 

We now have a two-parameter (rib, N) family of models 7i for a fixed realization of the 
Jij. It is very instructive to consider the phase diagram of H as a function of these two 
parameters (Fig. 0). 

The limit of large rib, with N fixed leads to the semi-classical theory. For the special case of 
SU(2) antiferromagnets with a two-sublattice collinear Neel ground state, the semiclassical 
fluctuations are described by the 0(3) non-linear sigma model of Sections 3,4. For other 
models \\ 



BQ, 52, 53, p4l M, M, the structure of the non-linear sigma models is rather 



more complicated and will not be considered here. 
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A second limit in which the problem simplifies is N large at fixed rib [571, [26| . It can be 



shown that in this limit the ground state is quantum disordered. Further, the low-energy 



dynamics of Tt is described by an effective quantum-dimer model [5S, g6|, with each dimer 
configuration representing a particular pairing of the sites into valence-bonds. This model 
is itself described by a non-trivial many-body Hamiltonian and has so far only be studied 



for the case of the unfrustrated square lattice antiferromagnet |p9| . p5| . |60| . |61f . We will not 
discuss it any further here, but existing results will be briefly noted later. 

The most interesting solvable limit is obtained by fixing the ratio of n& an d N 

k = £ (108) 



and subsequently taking the limit of large N [51]; this limit will be studied in this section 



in considerable detail. The implementation of 7i in terms of bosonic operators also turns 
out to be naturally suited for studying this limit. The parameter k is arbitrary; tuning k 
modifies the slope of the line in Fig. |7] along which the large N limit is taken. From the 
previous limits discussed above, one might expect that the ground state of 7i has magnetic 
LRO for large k and is quantum-disordered for small k. We will indeed find below that for 
any set of there is a critical value of k = k c which separates the magnetically ordered 
and the quantum disordered phase. 

A powerful feature of the bosonic large- N limit noted is the existence of a second-order 
phase transition at N = oo. In the vicinity of the phase transition, we expect the physics 
to be controlled by long-wavelength, low-energy spin fluctuations; the large- N method offers 
an unbiased guide in identifying the proper low-energy degress of freedom and determines 
the effective action controlling them. Having obtained a long-wavelength continuum theory 
near the transition, one might hope to analyze the continuum theory independently of the 
large- iV approximation and obtain results that are more generally valid. 

We will discuss the structure of the N = oo mean-field theory , with n& = kN in Section 
6. A. The long-wavelength effective actions will be derived and analyzed in Section 6.B. 
Finally topological Berry phase effects will be considered in Section 6.C. 

6. A Mean-field theory 

We begin by analyzing TC at N = oo with rib — k>N. As noted above, this limit is most 
conveniently taken using the bosonic operators. We may represent the partition function of 
H by 

Z = [ VQVbVXexp I- [ Cdr) , (109) 
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where 
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J -^J aP b^ + H.c. 



(110) 



Here the Aj fix the boson number of rib at each site; r-dependence of all fields is implicit; Q 
was introduced by a Hubbard-Stratonovich decoupling of TC. An important feature of the 
lagrangian £ is its U(l) gauge invariance under which 

bt 



Qi,j 



&L(0 ex P {ipi{r)) 
Q itj exp (-zpi(r) - ipj(r)) 
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The functional integral over C faithfully represents the partition function as long as we fix a 
gauge, e.g. by the condition dX/dr — at all sites. 

The 1/N expansion of the free energy can be obtained by integrating out of L the 2iV- 
component b,b fields to leave an effective action for Q, X having co-efficient iV (since oc N). 
Thus the iV — ► oo limit is given by minimizing the effective action with respect to "mean- 
field" values of Q = Q, X = X (we are ignoring here the possibility of magnetic LRO which 
requires an additional condensate x a = (b a ) - this has been discussed elsewhere [§, 0, §). 
This is in turn equivalent to solving the mean-field Hamiltonian 



n 



MF 



E 

<i,j> 



N 
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J ijQ 
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This Hamiltonian is quadratic in the boson operators and all its eigenvalues can be deter- 
mined by a Bogouibov transformation. This leads in general to an expression of the form 



Hmf = E MF [Q, X]+Y. ^[0, MX 



(113) 



The index \i extends over 1 . . .number of sites in the system, Emf is the ground state energy 
and is a functional of Q, A, uo^ is the eigenspectrum of excitation energies which is a also a 
function of Q, A, and the 7^ represent the bosonic eigenoperators. The excitation spectrum 
thus consists of non-interacting spinor bosons. The ground state is determined by minimizing 
Emf with respect to the Qij subject to the constraints 

dE M F 



dX,. 







The saddle-point value of the Q satisfies 

Qij = (J a pK^) 



(114) 



(115) 
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Note that Qij = —Qji indicating that Qij is a directed field - an orientation has to be chosen 
on every link. 

We will now consider the ground state configurations of the Q, A fields and the nature of 
the bosonic eigenspectrum for a variety of non-random antiferromagnets: 

Q.A.I Ji — J 2 — J3 model 

This is the square lattice antiferromagnet with first (Ji), second (J 2 ), and third ( J3) neighbor 



interactions [52], [53j]. We examined the values of the energy Emf for Qij configurations 
which had a translational symmetry with two sites per unit cell. For all parameter values 
configurations with a single site per unit cell were always found to be the global minima. 
We will therefore restrict our attention to such configurations. The \ field is therefore 
independent of i, while there are six independent values of Q^: 

Qi,i+x = Ql,x 

Qi,i+y = Ql,y 

Qi,i+y+x = Ql,y+x 

Qi,i+y—x = Ql,y—x 

Qi,i+2x = Q3,x 

Qi,i+2y = Qz,y (H6) 

For this choice, the bosonic eigenstates are also eigenstates of momentum with momenta k 
extending over the entire first Brillouin zone. The bosonic eigenenergies are given by 



|2 



1/2 



CJ k = (A 2 - \Ay 

A k = J\ (Q ljX sin k x + Q 1)V sin k y ) + J 2 (Q 2 ,y+ X sm(k y + k x ) + Q 2 , y - X sin(/c y - k x )) 

+ J3 (Qs,x sin(2^) + Q 3tV sm(2k y )) (117) 

We have numerically examined the global minima of Emf as a function of the three 
parameters J 2 / Ji, J3/J1, and N/rib 0, 0|. The values of the Qij at any point in the phase 
diagram can then be used to classify the distinct classes of states. The results are summarized 
in Figs. ||-11 which show various sections of the three-dimensional phase diagram. All of the 



phases are labeled by the wavevector at which the spin structure factor has a maximum. 
This maximum is a delta function for the phases with magnetic LRO, while it is simply a 
smooth function of k for the quantum disordered phases (denoted by SRO in Figs j3 HTT|) . The 
location of this maximum will simply be twice the wavevector at which cuk has a mimimum: 
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Figure 8: Ground states of the J\ — J2 — J3 model for J 3 = as a function of J2/J1 and 
N/n b (n b — 25* for SU(2)). Thick (thin) lines denote first (second) order transitions at 
N = 00. Phases are identified by the wavevectors at which they have magnetic long-range- 
order (LRO) or short-range-order (SRO). The links with Q p 7^ in each SRO phase are 
shown. The large N/n b , large J2/J1 phase has the two sublattices decoupled at N = 00; 
each sublattice has Neel -type SRO. Spin-Peierls order at finite N for odd n b is illustrated by 
the thick, thin and dotted lines. The (n,7i)-SR0 and the "decoupled" states have line-type 
spin-Peierls order for n b = 2 (mod 4) and are valence-bond-solids for n b = O(mod 4). The 
(0, tt)-SRO state is a valence-bond-solid for all even n b . 
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Figure 9: As in the previous figure but as a function of J 2 /Ji and J3/J1 for (a) N/rib = 1 
and (b ) N/n-b = 5. The inset in (a ) shows the region at the tip of the arrow magnified by 20: 
a direct first-order transition from {n,-n)-LRO to (0,tt)-LRO occurs up to J3/J1 = 0.005. 
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Figure 10: As in the previous figure but for J3/J1 = 0.35 
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Figure 11: As in the previous figure but for J3/J2 = 0.5 
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this is because the structure factor involves the product of two bosonic correlation functions, 
each of which consists of a propagator with energy denominator u^. 

Each of the phases described below has magnetic LRO for large n^/N and is quantum 
disordered for small n^/N. The mean-field result for the structure of all of the quantum 
disordered phases is also quite simple: they are featureless spin fluids with free spin- 1/2 
bosonic excitations ("spinons") with energy dispersion which is gapped over the entire 
Brillouin zone. Notice that this result is quite different from that of the O(M) non-linear 
sigma model of Sections 3,4 which found only integer spin excitations - the difference will be 
reconciled later. Some of the quantum disordered phases possess a broken lattice rotation 
symmetry even at iV = oo - these will be described below. The mininum energy spinons lie 
at a wavevector ko and decreases as n^/N. The onset of magnetic LRO occurs at the 
value of rib/N at which the gap first vanishes: co>k = 0. At still larger values of rib/N, we get 
macroscopic bose condensation of the b quanta at the wavevector ko, leading to magnetic 
LRO at the wavevector 2ko. 

We now turn to a description of the various phases obtained. They can be broadly 
classified into two types: 

Commensurate collinear phases 

In these states the wavevector ko remains pinned at a commensurate point in the Brillouin 
zone, which is independent of the values of J 2 1 Ji, J3/ J\ and n^/N . In the LRO phase the spin 
condensates on the sites are either parallel or anti-parallel to each other, which we identify 
as collinear ordering. This implies that the LRO phase remains invariant under rotations 
about about the condensate axis and the rotation symmetry is not completely broken. 
Three distinct realizations of such states were found 

1. (7T,7r) 

This is the usual two-sublattice Neel state of the unfrustrated square lattice and its quantum- 
disordered partner. These states have 

Ql,x — Ql,y 7^ 0) Q2,y+x = Ql,y-x = Q?,,x = Qz,y = (H8) 

From ( |117[ ), the minimum spinon excitation occurs at ko = ±(7r/2, 7r/2). The SRO states 
have no broken symmetry at N = oo. The boundary between the LRO and SRO phases 
occurs at N/rib < 2.5, independent of J2/J1 (Fig [8]). This last feature is surely an artifact of 
the large N limit. Finite N fluctuations should be stronger as J2/J1 increases, causing the 
boundary to bend a little downwards to the right. 
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2. (?r,0) or (0,tt) 
The (0, 7r) states have 



= 0, Q x ,y ^ 0, Q 2 , y +x = Q% y -x 7^ 0, and Q 3)X = Q 3>2/ = (119) 

and minimum energy spinons at k = ±(0,7r/2). The degenerate (vr, 0) state is obtained 
with the mapping x <-> y. The SRO state has a two-field degeneracy due to the broken 
x <-» y lattice symmetry. The LRO state again has two-sublattice collinear Neel order, but 
the assignment of the sublattices is different from the (ir, tt) state. The spins are parallel 
along the x-axis, but anti-parallel along the y-axis. 

An interesting feature of the LRO state here is the occurence of order- from-disorder [52 1 . 



Recall that the classical limit {n^/N = oo) of this model has an accidental degeneracy 
for J2/J1 > 1/2: the ground state has independent collinear Neel order on each of the A 
and B sublattices, with the energy independent of the angle between the spins on the two 
sublattices. Quantum fluctuations are included self-consistently in the N = 00, n b /N finite, 
mean- field theory, and lead to an alignment of the spins on the sublattices and LRO at (0, tt). 
The orientation of the ground state has thus been selected by the quantum fluctuations. 

The (0, 7r) states are separated from the (tt, tt) states by a first-order transition. In 
particular, the spin stiffnesses of both states remain finite at the boundary between them. 
This should be distinguished from the classical limit in which the stiffness of both states 
vanish at their boundary J 2 = Ji/2; the finite spin stiffnesses are thus another manifestation 
of or der-from- disorder. At a point well away from the singular point J 2 = Ji/2, n^/N = 00 
in Fig ^ the stiffness of both states is of order N(rib/N) 2 for N = 00 and large n^/N; near 
this singular point however the stiffness is of order Nfob/N) is induced purely by quantum 
fluctuations. These results have since also been obtained by a careful resummation of the 
semiclassical expansion IBBI, 153 . 



3. "Decoupled" 

For J2/J1 and N/n^ both large, we have a "decoupled" state (Fig H) with 

Q2, y+X = Q2, v -x ^ and Q 1 = Q 3 = 0. (120) 

In this case Q p is non-zero only between sites on the same sublattice. The two sublattices 
have Neel type SRO which will be coupled by finite N fluctuations. The N = 00 state does 
not break any lattice symmetry. This state has no LRO partner. 
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Incommensurate phases 

In these phases the wavevector ko and the location of the maximum in the structure factor 
move continuously with the parameters j65| . The spin-condensate rotates with a period 
which is not commensurate with the underlying lattice spacing. Further the spin condensate 
is coplanar: the spins rotate within a given plane in spin space and are not collinear. There 
is this no spin rotation axis about which the LRO state remains invariant. 

Further, no states in which the spin condensate was fully three dimensional ("double- 
spiral" or chiral states) were found; these would be associated with complex values of Q p . All 
the saddle points possesed a gauge in which all the Q p were real. Time-reversal symmetry 



was therefore always preserved in all the SRO phases of Figs [8|-|Il . 

The incommensurate phases occur only in models with a finite J3 (Figs |9|-|TT|). There 
were two realizations: 

1. (ir,q) or (q,ir) 

Here q denotes a wavevector which varies continuously between and ir as the parameters 
are changed. The (q, tt) state has 

Qi,x 7^ Qi, y 7^ 0, Q 2 ,x+ y = Q2, y -x ^ 0, Q 3tX ^ and Q^ y = 0; (121) 

the degenerate (n, q) helix is obtained by the mapping x y. The SRO state has a two-fold 
degeneracy due to the broken x *-> y lattice symmetry. 

2. (q,q) or (q,-q) 
The (q, q) state has 

Ql,x = Ql,y 7^ 0) Q2,x+y 7^ 0) Q2,y-x = 0, Qz,x = Q%,y 7^ 0j (122) 

this is degenerate with the (q, —q) phase and SRO state therefore has a two-fold degeneracy 
due to a broken lattice reflection symmetry. 

Note that the broken discrete symmetries in states with SRO at (0, 7r) and (q, ir) are 
identical: both are two-fold degenerate due to a breaking of the x y symmetry. The states 
are only distinguished by a non-zero value of Q3 in the (q, tt) phase and the accompanying 
incommensurate correlations in the spin-spin correlation functions. However Q 3 is gauge- 
dependent and so somewhat unphysical as an order parameter. In the absence of any further 
fluctuation-driven lattice symmetry breaking (see below), the transition between SRO at 
(0, 7i") and (q, n) is an example of a disorder line [p6fl ; these are lines at which incommensurate 
correlations first turn on. 
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An interesting feature of Figs 10-11 is that the commensurate states squeeze out the 
incommensurate phases as N/rib increases in both phase diagrams. We expect that this sup- 
pression of incommensurate order by quantum fluctuations is a general feature of frustrated 
antiferromagnets. This result is also consistent with the natural hypothesis that the states of 
the large N, fixed, but small n^/N should be consistent with the states of the large N, fixed 
rib theory (Fig. ^). This latter limit is described by the quantum-dimer model [5S] which is 
necessarily associated only with commensurate states. 

6. A. 2 Triangular and Kagome Lattices 

We have also examined the mean-field equations for the nearest-neighbor antiferromagnet 
on the triangular |)7| and kagome lattices in considerable detail ||. In both cases we found 



that the SRO phases had the full symmetry of the underlying lattice. They differ in this 
manner from all of the SRO phases of the square lattice. Further the Qij were real on every 
link and had the same magnitude \Q%j\ = Q- The only remaining degree of freedom is that 
associated with assigning an orientation to each link: minimization of the energy determined 
a unique orientation upto gauge-equivalent configurations. 

On the triangular lattice the spectrum of the free spin-1/2 bosonic spinon excitations 
was found to be given by 

1 /2 

o;(k) = (A 2 - J 2 Q 2 (sin ki + sin k 2 + sin fc 3 ) 2 ) (123) 

where the momentum k ranges over the first Brillouin zone of the triangular lattice and 

h p = k • e p (124) 

with the e p being unit vectors of length a pointing along the directions of the bonds on the 
triangular lattice 

ei = a(l/2,v / 3/2) 
e 2 = a(l/2, -V3/2) 

e 3 = a(-l,0) (125) 

A surface plot of this spectrum is shown in Fig [12] for nb/N = 0.25. The minimum-energy 
spinons are those at k i = (47r/3a)(l, 0) and k 2 = (47r/3a)(— 1, 0) and other points separated 
from these points by vectors of the reciprocal lattice generated by Gi = (47r/v / 3a)(0, 1) and 
G 2 = (47r/v / 3a)(v / 3/2, —1/2). For n b /N > 0.34 this state acquires magnetic LRO with the 
conventional three-sublattice Neel ordering of the triangular lattice. The spins are coplanar 
pointing towards the vertices of an equilateral triangle. 
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Figure 12: Momentum dependence of the energy, u(k) of the lowest excited spinon state 
for the quantum-disordered ground state of the triangular lattice quantum antiferromagnet 
at k — 0.25. The minimum excitations are the spinons at k = k x = (47r/3a)(l, 0) and 
k = k 2 = (47r/3a)(— 1, 0) and other points separated from k 1? k 2 by vectors of the reciprocal 
lattice 
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Very similar results were obtained on the kagome lattice. As this is not a Bravais lattice, 
it is necessary to introduce three sites per unit cell although no symmetry of the lattice is 
broken by the SRO state. To determine the free spin-1/2 bosonic excitation spectrum we 
need the following matrix 

(0 sin&i sin&;3 ^ 
sinfci sinA^ . (126) 
sin &3 sin k-z / 

(the k p were defined in Eqs. ( |1 24| ) and (|125| )) and to solve the eigenvalue equation 

pt(k)P(k)^(k)=pJ(k)^(k) (127) 

where p 2 , \i = 1,2,3 are the eigenvalues and </? M (k) the eigenvectors. Finally the bosonic 
eigenspectrum is given by 

^(k) = (A 2 -pJ(k)) 1/2 (128) 

The lowest energy spinon spectrum for n^/N = 0.35 is plotted in Fig [13]. The minimum 
energy spinons are at k = ki = (27r/3a)(l, 0) and k = k 2 = (27r/3a)(— 1, 0) and other points 
separated from k 1; k 2 by vectors of the reciprocal lattice generated by Gi = (2ir/ y / 3a)(0, 1) 
and G 2 = (27r/\/3a)(v^3/2, — 1/2). There is no quite good evidence that the spin-1/2, 



nearest-neighbor, kagome antiferromagnet is quantum disordered |68|, [69(]. The present cal- 
culation should therefore be relevant for this system. 

The spinon gap vanishes and magnetic LRO appears for nb/N > 0.53. The magnetic LRO 
is the configuration identified as the coplanar "v^3 x V3" structure in the literature |70, 71 



|72| , |73| . The huge accidental degeneracy of the classical ground state on the kagome lattice 
is completely lifted by the quantum fluctuations, and unique (upto global spin rotations) 
magnetic structure is obtained as the ground state. Note the rather natural way in which 
this happened directly in the the mean-field theory. The present large- iV approach thus seems 
to be ideally suited to examining order-from-disorder issues in frustrated antiferromagnets. 

6.B Fluctuations - long wavelength effective actions 

We now extend the analysis of Section 6.B beyond the mean-field theory and examine the 
consequences of corrections at finite N. The main question we hope to address are: 

• The mean-field theory yielded an excitation spectrum consisting of free spin-1/2 bo- 
sonic spinons. We now want to understand the nature of the forces between these 
spinons and whether they can lead to confinement of half-integer spin excitations. 

• Are there any collective excitations and does their dynamics modify in any way the 
nature of the ground state ? 
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Figure 13: Momentum dependence of the energy, u;(k) of the lowest excited spinon state 
for the quantum- disordered ground state (which has Q\ = —Q2) of the kagome lattice 
quantum antiferromagnet at k — 0.35. The energy is measures in units of J/2, and a is the 
nearest-neighbor spacing on the kagome lattice. The minimum excitations are the spinons 
at k = ki = (27r/3a)(l, 0) and k = k 2 = (27r/3a)(— 1, 0) and other points separated from 
ki , k 2 by vectors of the reciprocal lattice 
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The structure of the fluctuations will clearly be determined by the low-energy excitations 
about the mean-field state. We have already identified one set of such excitations: spinons at 
momenta near mimima in their dispersion spectrum, close to the onset of the magnetic LRO 
phase whence the spinon gap vanishes. An additional set of low-lying spinless excitations 
can arise from the fluctuations of the Qij and Aj fields about their mean-field values. The 
gauge-invariance ( |111|) will act as a powerful restriction on the allowed in the effective action 
for these spinless fields. We anticipate that the only such low-lying excitations are associated 
with the Aj and the phases of the Qij. We therefore parametrize 

Qi,i+e p = Qi,i+e p exp (-iQp) (129) 

where the vector e p connects the two sites of the lattice under consideration, Q is the mean- 
field value, and B p is a real phase. The gauge invariance (|1 1 1| implies that the effective 
action for the G p must be invariant under 

©p + Pi + Pi+e p - (130) 

Upon performing a Fourier transform, with the link variables P placed on the center of the 
links, the gauge invariance takes the form 

6 p (k) -> e p (k) + 2p(k) cos(£; p /2) (131) 

where k p = k-e p . This invariance implies that the effective action for the G p , after integrating 
out the b quanta, can only be a function of the following gauge-invariant combinations: 

I pg = 2 cos( V2)0 p (k) - 2 cos(k p /2)e q (k) (132) 

We now wish to take the continuum limit at points in the Brillouin zone where the action 
involves only gradients of the G p fields and thus has the possibility of gapless excitations. 
This involves expanding about points in the Brillouin zone where 

cos(fc p /2) = for the largest numbers of e p (133) 

We will apply this general principle to the models considered in Sec. 6. A. 

6.B.1 J x -J 2 - h Model 



We begin by examining the (7r,7r)-SRO phase. As noted in ( |118|) , this phase has the mean 



field values Qi tX = Q\ jV ^ 0, and all other zero. Thus we need only examine the condition 



( |133|) with e p = e x ,e y . This uniquely identifies the point k = G = (n, n) in the Brillouin 
zone. We therefore parametrize 



9 x (r) = ie iG - r A x (v) (134) 
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and similarly for 0^; it can be verified that both G and A x are real in the above equation. 
We will also be examining invariances of the theory under gauge transformations near G: so 
we write 



P(r) 



AG-r 



<p(r) 



(135) 



It is now straightforward to verify that the gauge transformations ( |131|) are equivalent to 

A x ^A x + d x ip (136) 

and similarly for A y . We will also need in the continuum limit the component of A near the 
wavevector G. We therefore write 



\ i = \ + ie %Gr A T {v l 



Under gauge transformations we have 



A T — > A T + d T tp 



(137) 



'138) 



Thus A x , A y , A T transform as components of a continuum U(l) vector gauge field. 

We will also need the properties of the boson operators under if. From ( 111 ) and (|135|) 
we see that the bosons on the two sublattices (A, B) with opposite charges ±1: 



b B 



b A ^ 
b B e- iip 



(139) 



Finally, we note that the bosonic eigenspectrum has a minimum near k = k = (7r/2,7r/2); 
we therefore parametrize 



°Ai 



(140) 



We insert the continuum parametrizations ( 134 ), (|137|) and (140) into the functional 
integral (|110|) , perform a gradient expansion, and transform the Lagrangian £ into 



(A + Mr) W + 1>? (J: ~ Mr) ^ a + A (Kl 2 + |^ a | 2 ) 
- UiQ x + VkVT) + JiQia 2 [(V + iA) ^ (V - iA) & 



2a 



+ (V - iA) ri a (V + iA) 



(141) 
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We now introduce the fields 



z°> = «+^*)/v / 2 

From Eqn ( |141| ), it is clear that the the n fields turn out to have mass A + 4J\Qi, while the 
z fields have a mass A — 4 J\Q\ which vanishes at the transition to the LRO phase. The tt 
fields can therefore be safely integrated out, and L yields the following effective action, valid 
at distances much larger than the lattice spacing p4|, |35|: 



S ^-IWaC if {^-'' A ^ f + ^% <142) 

Here /i extends over x,y,z, c = \^8JiQia is the spin-wave velocity, A = (A 2 — WJ^Qf) 1 ^ 2 
is the gap towards spinon excitations, and Af = A T /c. Thus, in its final form, the long- 
wavelength theory consists of a massive, spin-1/2, relativistic, boson z a (spinon) coupled to 
a compact U(l) gauge field. 

At distances larger than c/A, we may safely integrate out the massive z quanta and 
obtain a a compact U(l) gauge theory in 2+1 dimensions. This theory was argued by 



Polyakov [36, 74] to be permanently in a confining phase, with the confinement driven by 
"instanton" tunnelling events. The compact U(l) gauge force will therefore confine the 
z a quanta in pairs. In the present theory, the confinementlength scale turns out to be 
exponentially large in N: ~ e cN where the constant c diverges logarithmically as A — > [46|. 
There are thus no free spin-1/2 bosonic excitations for any finite N and all low-lying modes 
carry integral spin. The presence of an unbroken U(l) gauge force in the fluctuations has 
therefore completely disrupted the simple mean-field structure of these states. The spectrum 
of the quanta also acquires a gap from instantons effects which is exponentially small in 



N [T741 IB] 



The properties of the (0, 7r) phase are very similar to those of the (tt, tt) phase considered 
above, and will therefore not be discussed here. It can be shown quite generally that any 
quantum disordered state which has appreciable commensurate, collinear spin correlations 
willl have similar properties: confined spinons and a collective mode described by a compact 
£7(1) gauge theory. 

We now turn to a study of the incommensurate phases. It is not difficult to show that 



in this case it is not possible to satisfy the constraints ( 133 ) at any point in the Brillouin 
zone for all the non-zero Q p . This implies that there is no gapless collective mode in the 
incommensurate SRO phases. The structure of the theory is simplest in the vicinity of a 
transition to a commensurate collinear phase: we now examine the effective action as one 



46 



moves from the (tt, 7r)-SRO phase into the (q, g)-SRO phase (Figs 10-11]) (a very similar anal- 
ysis can be performed at the boundary between the (tt, 7r)-SRO and the (tt, q)- SRO phases). 
This transition is characterized by a continuous turning on of non-zero values of Qi^ + y + x, 



Qi,i+2x and Qi,i+2y- It is easy to see from Eqn ( |1 1 1|) that these fields transform as scalars 
of charge ±2 under the gauge transformation associated with A^. Performing a gradient 
expansion upon the bosonic fields coupled to these scalars we find that the Lagrangian £ of 
the (tt, 7r)-SRO phase gets modified to 

C^C+ f—($ A - (j a ^V^) + $ B • (j a ^2a^2p) + C.C.) (143) 

where $a,b are two-component scalars = (JzQz,x + J2Q2M+X, JzQz,y + J2Q2,y+x) with the 
sites on the ends of the link variables on sublattices A, B. Finally, as before, we transform 
to the z, tt variables, integrate out the tt fluctuations and obtain 

S eff = J J* df{\(8, - iA^z a \ 2 + r\z a \ 2 + I • (jafnTVzP) + c.c. + V($)} + . . . , 

(144) 

Here r = A 2 /c 2 , $ = ($^ + <&* B )/(2JiQia) is a scalar of charge —2; terms higher order 
in $ have been dropped. We have also added a phenomenological potential V($) which is 
generated by short wavelength fluctuations of the b a quanta. This effective action is also the 
simplest theory that can be written down which couples a spin- 1/2, charge 1, boson z a , a 
compact U(l) gauge field A^, and a two spatial component, charge —2, spinless boson $. It 
is the main result of this section and summarizes essentially all of the physics we are trying 
to describe. We now describe the various phases of S e ff 

1. Commensurate, collinear, LRO: (z a ) ^ 0, ($) = 

This is the state with commensurate, collinear, magnetic LRO 

2. Commensurate, collinear, SRO: (z a ) = 0, ($) = 

This is the quantum-disordered state with collinear spin correlations peaked at (tt, tt). 
Its properties where described at length above. The compact U(l) gauge force confines 
the z a quanta. The spinless collective mode associated with the gauge fluctuations also 
has a gap. 

3. Incommensurate, coplanar, SRO: (z a ) = 0, ($) ^ 

This is the incommensurate phase with SRO at (q, q) which we want to study. It is easy 
to see that condensation of $ necessarily implies the appearance of incommensurate 
SRO: ignore fluctuations of $ about ($) and diagonalize the quadratic form controlling 
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the z a fluctuations; the minimum of the dispersion of the z a quanta is at a non-zero 
wavevector 

k = (%))/2 (145) 

The spin structure factor will therefore have a maximum at an incommensurate wave- 
vector. This phase also has a broken lattice rotation symmetry due to the choice of 
orientation in the x — y plane made by $ condensate. The condensation of $ also has 
a dramatic impact on the nature of the force between the massive z a quanta. Detailed 
arguments have been presented by Fradkin and Shenker [75] that the condensation 



of a doubly charged Higgs scalar quenches the confining compact £7(1) gauge force 
in 2+1 dimensions between singly charged particles. Applied to the present problem, 
this implies that the charge —2 field $ condenses and deconfines the z a quanta. The 
excitation structure is therefore very similar to that of the mean-field theory: spin- 1/2, 
massive bosonic spinons and spinless collective modes which have a gap. The collective 
mode gap is present in this case even at N = oo and is associated with the condensation 
of $. 

4. Incommensurate, coplanar, LRO: (z a ) ^ 0, ($) ^ 

The condensation of the z quanta at the wavevector k above leads to incommensurate 
LRO, with the spin condensate spiraling in the plane. 

6.B.2 Triangular and Kagome Lattices 

We will focus mainly on the nature of the fluctuations on the triangular lattice M. The 
properties of the kagome lattice are very similar and have been discussed elsewhere ||. 

The magnetically ordered state on the triangular lattice is coplanar. From the results 
on the J\ — J2 — J3 model we anticipate that the fluctuations about the quantum disodered 
state on the triangular lattice will be similar to those of the incommensurate, coplanar, SRO 
states - there will be no gapless gauge modes and the spin-1/2 spinons remain unconfined. 

We now present a few details verifying this conjecture. It is not difficult to see that 



the constraint ( |133| ) for low-lying gauge modes can be satisfied at most two of the values of 
p — 1, 2, 3 (corresponding to the e p in Eqns. ( |125| )) at any point of the Brillouin zone. One 
such point is the wavevector 

2vr . . 
(0, 1) (146) 



where 

9al = 7T 
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9a3 = (147) 

Taking the continuum limit with the fields varying with momenta with close to g a we find 
that the I pq in ( |132|) depend only upon gradients of 0j and G>2- Under gauge transformations 
near the momentum g a , the bosons 

bf carry charges exp(zg a • r^). (148) 

It can be verified that these charges only take the values ±1 on the lattice sites. We have 
therefore imposed a certain 'staggering' of the charge assignments of the bosons which is 
quite analogous to that in the square-lattice antiferromagnets in Sec. 6.B.I. It is also helpful 
to parametrize the P in the following suggestive manner 

ex(r) = iA al {v)e^ 
e 2 (r) = -iA a2 (v)e^ r 

e 3 (r) = $ a (r)e lg - r (149) 

It can be verified that the condition for the reality of G p is equivalent to demanding that 
Aal, A a2 , $ a be real. We will now take the continuum limit with Aal, A a2 , $ a varying slowly 
on the scale of the lattice spacing. It is then not difficult to show that the invariants I pq then 
reduce to (after a Fourier transformation): 

In = e 2 • VA aX - ei • V ' A a2 

1 31 = e l -V$ a -2A al 

1 32 = e 2 -V§ a -2A a2 (150) 

Thus the A al , A a2 are the components of the connection of a gauge symmetry denoted U a (l); 
the components are taken along an 'oblique' co-ordinate system defined by the axes ii,e 2 . 
The field $ a transforms as the phase of charge 2 Higgs field under U a (l). 

A very similar analysis can be carried out near the two other points in the Brillouin zone 
where the other pairs of values of cos(fc p /2) vanish. These are the points 

2tt fy/3 1\ 



(151) 




bc x/3aV 2 ' 

which introduce the continuum symmetries £4(1) and U c (l) respectively. The P now reduce 
in the continuum limit to fields Ab 2 , A^ and A c i,§ c , A c3 respectively. Thus in the con- 
tinuum limit the lattice U(l) gauge symmetry has been replaced by a U a (l) x £4(1) x U c (l) 
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gauge symmetry. The three gauge symmetries correspond to the three different ways the 
triangular lattice can be distorted into a rectangular lattice with diagonal bonds: the phases 
on the horizontal and vertical bonds behave like gauge connections while the phases on 
the diagonal bonds become charge 2 Higgs fields. The system also possesses spin-1/2 bose 
excitations which carry charges ±1 of all 3 symmetries. 

The condensation of all of the Higgs fields is implicit, and there are therefore no low- lying 
physical gauge excitations. As in Sec. 6.B.1, we conclude that the instantons are quenched 
and that unit charges are expected to be unconfined ||75|| ; in particular the spin-1/2 bose 
quanta, which carry the U a (l) charges specified in ( |148| ), and analogous £4(1) and U C {1) 
charges, will remain unconfined. 

6.C Berry Phases 

We found two different radically different consequences of gauge fluctuations in the quantum- 
disordered phases in section 6.B: 

(a) the commensurate collinear phases appeared to have gapless gauge modes; however, they 
acquired a gap from instanton effects which also led to the confinement of spinons. 

(b) the coplanar phases on the square, triangular, and kagome phases had only gapped 
gauge fluctuations; the mean-field spectrum was stable to fluctuation effects and the bosonic 
spinons remained unconfined. 

It was also clear from the analyses in Sections 6.B and 3 that the Berry phases of small 
fluctuations about the mean-field saddle points had been properly accounted for. 

The only possibility that has not been explored carefuly is that the topologically non- 
trivial gauge field configurations (instantons, vortices) might posses some non-trivial Berry 
phases. Such phases are of course not present in a theory with a simple Maxwell action 
for the gauge-field fluctuations; in the present theory the effective action is much more 
complicated and the bosons propogate in a complicated gauge-field background can acquire 
Berry phases. Such effects can, in principle, be present in both the collinear and coplanar 
SRO phases. However, in the vicinity of the transition to the LRO phases, instanton effects 
are not expected to be dominant in the coplanar states; in contrast the collinear phases have 
already been shown to be radically modified by instanton effects. 

Detailed calculations of instanton Berry phase effects have been performed in the collinear 



phases. These calculations are rather involved and the reader is referred to Ref. |35[ for 
further details: we will simply present the results here. It was first shown by Haldane |[33|| . 
using the non-linear sigma model formulation of Section 3 for fluctuations in the (tt, it) LRO 
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Figure 14: The A, B sublattices of the lattice of spins and the sublattices W, X, Y, Z of the 
dual lattice. 

state, that 'hedgehog' n configurations had the following Berry phase term in its action 

S B = i 7 ^- m sC ( 152 ) 

^ a 

where n& = 2S, m s is the charge of an hedgehog centered on the square lattice plaquette 
numbered s, and ( s takes the values 0, 1, 2, 3 on the W, X, Y, Z plaquette sublattices (Fig [14]); 
it can be shown that no symmetry is broken by this choice of phases, and that all observable 
correlation functions are invariant under the full lattice symmetry group (in the absence 
of any dynamical symmetry breaking, of course). Note that this phase is always unity for 
even-integer spins. 



Subsequently, Refs |34], [35[] examined these Berry phases directly in the commensurate, 
collinear SRO phases. It was argued that the remnants of the hedgehogs in the LRO phases 
were precisely the instantons in the compact U(l) gauge field of Section 6.B. The Berry phase 
of bosons propagating in an instanton gauge-field background was evaluated and found to 
be identical to (|152|) in the (tt, tt) SRO phase; m s is now the total flux emannating from the 
instanton divided by 2n. The result in the (0,7r) SRO phase was different [fj]]: 

S B = inn b Y ™>s [Rax] (153) 

s 

where [R sx ] is the integer part of R x , the x-coordinate of the plaquette s; this phase is unity 
for all integer spins. The result in the decoupled state was imply two copies of the (n, tt) 
SRO phase. 

An analysis of the dynamics of the instantons was then carried out B3, |?6[ . In the SRO 



phase the instantons interact with a Coulombic 1/R potential; the instanton plasma can 
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rib = 1,3 (mod 4) rib = 2 (mod 4) rib = (mod 4) 
(a) (b) (c) 

Figure 15: Symmetry of non-Neel ground states of H as a function of n&(mod4) with the 
minimum possible degeneracies of 4,2,1 respectively (n b = 2S for SU(2)). The full-dotted- 
blank lines represent different values of (Sj • Sj) on the links. 



therefore be mapped onto a dual sine-Gordon model in which the instanton Berry phases 
appear as frustrating phase-shifts in the arguments of the cosine term. Finally analysis of 
this sine-Gordon model showed that non-unity Berry phases led to spontaneous breaking of 
a lattice rotational symmetry through the appearance of spin-Peierls order. The nature of 
the spin-Peierls ordering depended strongly on the value of n;,(mod4). 



For the case of (n, n) SRO the spin-Peierls ordering is shown in Fig [15|. The valence- 
bonds align in columns (4-fold degenerate) for half-integer spins, and along lines (2-fold 
degenerate) for integer spins. Only for even integer spins is there no breaking of symmetry 



and we obtain a valence-bond-solid state [50]. It is interesting that an analysis of the 



quantum-dimer model for n& = 1 also found spin-Peierls order of the columnar type. 
Thus the N — > oo, rib = 1 theory agrees with the iV — > oo, nb/N fixed, n& = 1 (mod 4) 
theory. 

In the "decoupled" SRO phase, the above analysis applies to each sublattice separately, 
giving for rib — 1 (mod 4) the type of spin-Peierls correlations shown in Fig [8]. There is a 
total of 4 x 4 = 16 states for this case but coupling between the sublattices will reduce this 
to 8 states, all of one of the two types shown. The state with the 'dimers' parallel to one 
another has more possibilities for resonance using the J\ bonds and is likely to be the ground 
state. For nb = 2 (mod 4) , there will be 2x2/2 = 2 states, and for nb = (mod 4) , just 
one. 

For the (0, 7r) SRO state the spin-Peierls order of the type shown in Fig |8| for rib odd 
(half-integer spins), and a VBS state for n& even (integer spins). Combined with the choice 
(0, 7i") or (7T, 0) this gives degeneracies 2, 4, 2, 4 for rib = 0, 1, 2, 3 (mod 4) . 
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6.D Summary 



The above analysis has been rather involved, but the essential results are rather simple. Let 
us recall here the main properties of the quantum disordered phases found in two dimensional 
ant if erromagnet s . 

• Commensurate, collinear phases 

The spin structure factor has a well-defined maximum at a commensurate point in the 
Brillouin zone. The spinon excitations are massive and confined in pairs - there are no 
low-lying excitations with half-integral spins. There is a spinless collective mode which 
has a gap induced by instanton effects. Spin-Peierls order is present and its nature 
depends on the value of nfc(mod4) (for the square lattice). 

• Incommensurate, coplanar phases 

The spin structure factor has well-defined maxima at incommensurate points in the 
Brillouin zone. The spin-1/2 spinon excitations are massive, deconfined and carry 
bosonic statistics. The spinless collective mode has a gap induced by the condensation 
of a Higgs field: the magnitude of the condensate is proportional to the incommensu- 
ration. There is a broken lattice reflection symmetry associated with the choice of an 
axis about which the fluctuating spiral order is present. 

• Commensurate, coplanar fluids on the triangular and kagome lattices 

These phases from the incommensurate coplanar states only in that the peak of the 
structure factor is at a commensurate point, and there is no broken symmetry. These 
states thus appear to violate Laughlin's fractional quantization principle JT7| that all 
spin-1/2 excitations of featureless spin-fluids should posses fractional statistics: the 
spinons in the present theory are bosonic. 

6. Comparison with numerical and series results 



Many numerical [78] and series analyses JfJj have appeared on the J\ — J2 — J3 model with 



J3 = for the spin-1/2 £77(2) model i.e. N — 1, — 1 in the notation of this paper. They 
find (7r,7r)-LRO at small J2/J1, (vr, 0)-LRO at large J2/J1 and an intermediate SRO phase 
around J2/J1 = 1/2. This is in agreement (see Fig [8]) with our prediction in Section 6.B that 
the phase boundary between (it, 7r)-LRO and (tt, it)- SRO bends downwards at finite N with 



increasing J%j J\. Analyses of this intermediate phase at J2/J1 [[78], [Fi| shows clear evidence 
of columnar spin-Peierls ordering j34j, also in agreement with the results of Section 6.C. An 
additional intermediate phase with (0, 7r)-SRO has not been ruled out. 
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There has also been intensive work recently on the ground state of the SU(2) antifer- 
romagnet on the triangular lattice. For spin-1/2, there is good evidence |68], B9| that the 



ground state is quantum disordered. There are preliminary indications [80] that spin-Peierls 
order is absent in this state (which was predicted in Ref. and discussed in Section 6.B.2), 
although larger system sizes are required before any firm conclusion can be reached. Another 
interesting issue of the kagome is the nature of the large-spin magnetically ordered ground 



state selected by quantum fluctuations. A number of investigators |70|, [71], [T2|, [7IJ have found 
the \/3 x a/3 discussed in Section 6. A. 2. 

7. CONCLUSIONS 

It should be clear from all the discussion in this paper that there has been a great deal 
of progress in our understanding of non-random two dimensional antiferromagnets. The 
magnetically ordered states had been understood many years ago by semiclassical analyses. 
This course therefore focussed on recent work on the nature of the quantum phase transition 
to the quantum disordered phase and the properties of the quantum disordered phase itself. A 
partial classification of the different types of possible quantum disordered states has emerged, 
and was summarized in Section 6.D - the reader is urged to review this section, even if he/she 
did not have the patience to read the rest of Section 6. Also interesting was the fundamental 
connection between the properties of the quantum disordered states and the nature of the 
spin ordering in the magnetically ordered states. 

Comparison of these results with experiments is complicated by the ubiquitous pres- 
ence of randomness. Unlike classical systems, even weak randomness has a strong effect on 
the properties of quantum phase transitions and the low-energy properties of the quantum- 
disordered phase. In the renormalization group language, randomness is almost always a 
relevant perturbation at fixed points controlling the properties of quantum phase transitions 
and quantum-disordered phase. Our comparision with experiments has therefore been lim- 
ited to a simple scaling analysis of phase transitions in random quantum antiferromagnets - 
this was summarized in Section 5. More detailed study of random quantum antiferromagnets 
is therefore clearly a high priority for future work. Recently, there has been some progress 
in solving random quantum magnets with infinite-range exchange interactions [24, 25|; a 
variety of unusual results were obtained, including the presence of gapless excitations in the 
quantum disordered phase. There have also been exact solutions of random quantum spin 
chains which are very instructive One hopes that these works are just the first stages 
of much further work on random quantum spin systems: the outlook for more interesting 
results is promising and the number of available experimental systems continues to grow. 



54 



Acknowledgements 

I thank N. Read and Jinwu Ye for several fruitful collaborations which led to the work 
reviewed here. This research was supported by NSF Grant No. DMR-8857228 and by a 
fellowship from the A. P. Sloan Foundation. 

References 

[1] HA. Bethe, Z. Physik, 71 205 (1991). 

[2] P.W. Anderson, Phys. Rev. B 86, 694 (1952). 

[3] T. Nagamiya, R. Kubo, and K, Yosida, Phil. Mag., Suppl., 14, 1 (1955). 
[4] P.W. Anderson, Science 235, 1196 (1987). 

[5] J.G. Bednorz and KA. Muller, Z. Phys. B64, 188 (1986); M.K. Wu et. al, Phys. Rev. 
Lett. 58, 908 (1987). 

[6] N. Read and S. Sachdev, Phys. Rev. Lett. 66, 1773 (1991). 

[7] S. Sachdev and N. Read, Int. J. Mod. Phys. B5, 219 (1991). 

[8] S. Sachdev, Phys. Rev. B , 45, 12377 (1992). 

[9] S.M. Hayden, G. Aeppli, H. Mook, D. Rytz, M.F. Hundley, and Z. Fisk, Phys. Rev. 
Lett. , 66, 821 (1991). 

[10] S.-W. Cheong, G. Aeppli, T.E. Mason, H. Mook, S.M. Hayden, PC. Canfield, Z. Fisk, 
K.N. Clausen, and J.L. Martinez, Phys. Rev. Lett. , 67, 1791 (1991). 

[11] B. Keimer, R.J. Birgeneau, A. Cassanho, Y. Endoh, R.W. Erwin, MA. Kastner, and 
G. Shirane, Phys. Rev. Lett. , 67, 1930 (1991). 

[12] S.M. Hayden, G. Aeppli, R. Osborn, A.D. Taylor, T.G. Perring, S.-W. Cheong, and Z. 
Fisk, Phys. Rev. Lett. , 67, 3622 (1991). 

[13] T.E. Mason, G. Aeppli, and HA. Mook, Phys. Rev. Lett. , 68, 1414 (1992). 

[14] T.R. Thurston, R.J. Birgeneau, Y. Endoh, P.M. Gehring, MA. Kastner, H. Kojima, M. 
Matsuda, G. Shirane, I. Tanaka, and K. Yamada, MIT preprint. 



55 



[15] B. Keimer, N. Belk, R.J. Birgeneau, A. Cassanho, C.Y. Chen, M. Greven, M.A. Kastner, 
A. Aharony, Y. Endoh, R.W. Erwin, and G. Shirane, MIT preprint. 

[16] M. Matsuda, Y. Endoh, K. Yamada, H. Kojima, I. Tanaka, R.J. Birgeneau, M.A. Kast- 
ner, and G. Shirane, MIT preprint. 

[17] A. P. Ramirez, G.P. Espinosa, and A.S. Cooper, Phys. Rev. Lett. , 64, 2070 (1990). 

[18] C. Broholm, G. Aeppli, G.P. Espinosa, and A.S. Cooper, Phys. Rev. Lett. , 65, 3173 
(1990). 

[19] G. Aeppli, C. Broholm, and A. Ramirez in Proceedings of the Kagome Workshop, NEC 
Research Institute, Princeton NJ. 

[20] D.S. Greywall, Phys. Rev. B 41, 1842 (1990). 

[21] V. Elser, Phys. Rev. Lett. , 62, 2405 (1990). 

[22] S. Sachdev and Jinwu Ye, Phys. Rev. Lett. 69, 2411 (1992). 

[23] S. Chakravarty, B.I. Halperin, and D.R. Nelson, Phys. Rev. Lett. 60, 1057 (1988); Phys. 
Rev. B 39, 2344 (1989). 

[24] S. Sachdev and J. Ye, "Gapless spin-fluid ground state in a random quantum Heisenberg 
magnet", paper 9212030 on cond-mat@babbage.sissa.it. 

[25] J. Ye, S. Sachdev, and N. Read, "A solvable spin-glass of quantum rotors", paper 
9212027 on cond-mat@babbage.sissa.it. 

[26] N. Read and S. Sachdev, Nucl. Phys. B316, 609 (1989). 

[27] Coherent Statesby J.R. Klauder and B. Skagerstam, World Scientific, Singapore (1985); 
see also Generalized Coherent States and Their Applications, by A. Perelomov, Springer- 
Verlag, New York (1986). 

[28] R. P. Feynman, Phys. Rev. 84, 108 (1951). 

[29] F.D.M. Haldane, Phys. Lett. 93A, 464 (1983); Phys. Rev. Lett. 50, 1153 (1983); and 
J. Appl. Phys. 57, 3359 (1985). 

[30] I. Affleck, Nucl. Phys. B257, 397 (1985); Nucl. Phys. B265, 409 (1985). 

[31] I. Affleck and F.D.M. Haldane, Phys. Rev B36, 5291 (1987). 



56 



[32] X.G. Wen and A. Zee, Phys. Rev. Lett., 61, 1025 (1988); E. Fradkin and M. Stone, 
Phys. Rev. B38, 7215 (1988); T. Dombre and N. Read, Phys. Rev. B38, 7181 (1988); 
R. Shankar and S. Sachdev, unpublished. 

[33] F.D.M. Haldane, Phys. Rev. Lett. 61, 1029 (1988). 

[34] N. Read and S. Sachdev, Phys. Rev. Lett. 62, 1694 (1989). 

[35] N. Read and S. Sachdev Phys. Rev. B 42, 4568 (1990). 

[36] A.M. Polyakov, Gauge Fields and Strings, Harwood, New York (1987). 

[37] S. Sachdev and Jinwu Ye, unpublished. 

[38] S-k. Ma, Modern Theory of Critical Phenomena, Benjamin/Cummings, Reading (1976). 
[39] R. Abe, Prog. Thoer. Phys. 49, 1877 (1973). 
[40] B. Keimer, private communication. 

[41] R.J. Gooding and A. Mailhot, Phys. Rev. B 44, 11852 (1991); A. Aharony et. al, Phys. 
Rev. Lett. 60, 1330 (1988). 

[42] K. Binder and A.P Young, Rev. Mod. Phys. 58, 801 (1986). 

[43] G.A. Baker et.al. Phys. Rev. B 17, 1365 (1978). 

[44] J.T. Chayes et.al. Phys. Rev. Lett. 57, 2999 (1986). 

[45] D.S. Fisher, private communication. 

[46] G. Murthy and S. Sachdev, Nucl. Phys. B344, 557 (1990). 

[47] S.N. Dorogovstev, Phys. Lett. 76A, 169 (1980); D. Boyanovsky and J.L. Cardy, Phys. 
Rev. B 26, 154 (1982); I.D. Lawrie and V.V. Prudnikov, J. Phys. C, 17, 1655 (1984). 

[48] C.C. Wan, A.B. Harris, and J. Adler, J. Appl. Phys. 69, 5191 (1991). 

[49] R.N. Bhatt and PA. Lee, Phys. Rev. Lett. 48, 344 (1982). 

[50] I. Affleck, T. Kennedy, E.H. Lieb, and H. Tasaki, Phys. Rev. Lett. 59, 799 (1987); 
D.Arovas, A. Auerbach and F.D.M. Haldane, Phys. Rev. Lett. 60, 531 (1988). 

[51] D.P Arovas and A. Auerbach, Phys. Rev. B38, 316 (1988); Phys. Rev. Lett. 61, 617 
(1988). 

57 



[52] T. Einarsson and H. Johannesson, Phys. Rev. B43, 5867 (1991); T. Einarsson, P. 
Frojdh, and H. Johannesson, Phys. Rev. 45, 13121 (1992). 

[53] P. Chandra, P. Coleman, and A.I. Larkin, Phys. Rev. Lett. 64, 88 (1990); P. Chandra 
and P. Coleman, Int. J. Mod. Phys. B3, 1729 (1989). 

[54] B.I. Halperin and W.M. Saslow, Phys. Rev. B 16, 2154 (1977). 

[55] T. Dombre and N. Read, Phys. Rev. B 39, 6797 (1989). 

[56] P. Azaria, B. Delamotte, and T. Jolicoeur, Phys. Rev. Lett. 64, 3175 (1990). 

[57] I. Affleck and J. Brad Marston, Phys. Rev. B37, 3774 (1988). 

[58] D. Rokhsar and S. Kivelson, Phys. Rev. Lett. 61, 2376 (1988). 

[59] S. Sachdev, Phys. Rev. B 40, 5204 (1989). 

[60] E. Fradkin and S. Kivelson, Mod. Phys. Lett. B4, 225 (1990). 

[61] S. Sachdev and R. Jalabert, Mod. Phys. Lett. B4, 1043 (1990). 

[62] C. Henley, Phys. Rev. Lett. 62, 2056 (1989). 

[63] A. Chubukov, Phys. Rev. B 44, 392 (1991). 

[64] F. Mila, D. Poilbanc, and C. Bruder, Phys. Rev. B 43, 7891 (1991). 

[65] Similar incommensurate states with valence bonds connecting sites on the same sublat- 
tice are also found in a study of doped systems; B.I. Shraiman and E.D. Siggia, Phys. 
Rev. Lett. 62, 1564 (1989); C. Jayaprakash, H.R. Krishnamurthy and S. Sarkar, Phys. 
Rev. B 40, 2610 (1989); C.L.Kane et al, ibid 41,2653 (1990); K. Feinsberg, P. Hedegard 
and M. Brix Pederson, ibid 40, 850 (1989). 

[66] J. Stephenson, Can. J. Phys. 48, 2118, 1724 (1970). 

[67] D. Yoshioka and J. Miyazaki, J. Phys. Soc. Jpn, 60, 614 (1991); A. Kol and N. Read, 
unpublished 

[68] C. Zeng and V. Elser, Phys. Rev. B 42, 8436 (1990). 

[69] R.R.P. Singh and D.A. Huse, Phys. Rev. Lett. , 68, 1766 (1992). 

[70] A. Chubukov, Phys. Rev. Lett. , 69, 834 (1992). 



58 



[71] D.A. Huse and A.D. Rutenberg, Phys. Rev. B 45, 7536 (1992). 

[72] A.B. Harris, C. Kallin, and A.J. Berlinsky, Phys. Rev. B 45, 2889 (1992). 

[73] J.T. Chalker, P.C.W. Holdsworth, and E.F. Shender, Phys. Rev. Lett. , 68, 855 (1992). 

[74] A.M. Polyakov, Nucl. Phys. B120, 429 (1977); 

[75] E. Fradkin and S.H. Shenker, Phys. Rev. D19, 3682 (1979). 

[76] W. Zheng and S. Sachdev, Phys. Rev. B 40, 2704 (1989). 

[77] R.B. Laughlin, Science 242, 525 (1988). 

[78] E. Dagotto and A. Moreo, Phys. Rev. Lett. 63 2148 (1989); R.R.P. Singh and 
R. Narayan, Phys. Rev. Lett. 65, 1072 (1990); A. Moreo, E. Dagotto, T. Joliceour, 
and J. Riera, Phys. Rev. B 42, 6283 (1990); T. Ziman and H. Schulz, preprint 

[79] M.P Gelfand, R.R.P. Singh, and D.A. Huse, Phys. Rev. B 40, 10801 (1989); 
M.P Gelfand, Phys. Rev. B 42, 8206 (1990). 

[80] J.T. Chalker and J.F.G. Eastmond, Oxford preprint, OUTP-92-23S. 

[81] C.A. Doty and D.S. Fisher, Phys. Rev. B45, 2167 (1992); D.S. Fisher, Phys. Rev. Lett. 
69, 534 (1992); unpublished. 



59 



